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A hypothetical aelfconsistent picture of evolution of prestellar interstellar 
dust through a comet phase leades to predictions about the composition of the 
' circum-solar dust cloud • Scattering properties of thus resulting conglomerates 
with a "bird's-nest" type of structure are investigated using a micro-wave 
analogue technique. Approximate theoretical methods of general interest are 
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Introductio n 

in mum mmam nmmumm mmm ***** 

The interplanetary particles producing the zodiacal light by scattering of 
sunlight were for many years considered to be in a range of sizes not very 
widely ’different from the wavelength of visible radiation. Furthermore, most 

I 

of the calculations of the brightness and polarization were done for rpherical 
particles using Mi# theory (for references see Weinberg, 1967, and the review 
article by Weinberg and Sparrow, 1278). A simple preliminary theoretical 
approach showed that significant modifications in tho prediction were possible, 
if one used similar-sized particles but non-spherical rather than spherical 
ones (Greenberg, 1970). 

It 

Recent observations of the optical properties (Weinberg and Sparrow, 1978) 
and the size and composition of the interplanetary particles as determined from 
collection and impact studies (Fechtig, 1976, and Brownlee, 1978) have led to 
new concepts of the interplanetary particles. Large irregular particles in the 
10 pm to 100 pm size range are now believed to provide most of the zodiacal 
light (Giese et al. , 1978). To date, only two models have reasonably well 
reproduced the zodiacal-light brightness and degree of polarization as deduced 
by Dumont (19?6) from observations. It is interesting that the two models were 

t 

arrived at using quite different approaches. Fluffy interplanetary dust grains 
were studied by Giese et al. (1978) with connection to the type of particles 
collected in the Earth's upper atmosphere (Brownlee, 1978). Our approach is to 
apply an evolutionary model leading to the concept that we call a "bird's-nent" 
structure (Greenberg and Gustafson, 1980). 
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Chop tor 1 

THE PURPOSE OF THE INVESTIGATION 

.The hypothesis to be tested in this study is n unifying evolutionary modal 
relating interstellar duat to comets and finally to the zodiacal-light particles » 
We will start with a set of assumptions that satisfy the existing observations 
but that may not be unique, It will be shown in chapter 5 that the resulting 
predictions are consistent with observations that in the absence of this model 
were unrelated. We adopt the bi-modtil model for interstellar dust of Greenberg 
and Hong (1973), The comets are assumed to have aggregated directly from the 
dust during the formation of the Solar System. Finally, the zodiacal-light 
particles are modeled as cometary debris. 

Although there is an identifiable mechanism for internal heating by 
primordial radioactive species, which could make a significant change in the 
internal structure of a comet nucleus (Irvine et al,, 1980), our modeling is an 
attempt to work as closely as possible with constructs based on the hypothesis 
that we are seeing pieces from the surfaces of essentially primordial eometa. 

One justification for doing this is that there is as yet no way of knowing to 
what extent comets of various sizes may have undergone internal or external 
structural or chemical modifications by heating. There are uncertainties both in 
the amount of primordial radioactive species accreted in the primordial comets 
as well as in the details of the heating process (d’Hendecourt and Greenberg, 
unpublished calculations) . Another justification is that the albedo of comets 
(Delaemme and Rud, 1973) and the abundance of volatiles seem more consistent with 
at least the outer parts of a significant number of comets remaining primordial 
than being subject to high hooting. Assuming that comets have undergone little 
change in their properties since birth - and this seems to be an inherent 
assumption in the '‘icy" models of Whipple (1978) and Delsemme and Rud (1973) - 
then the cometary debris could essentially be derivates of clumps of large 
numbers of interstellar grains. 

The major constituents in the. bi-modal distribution of interstellar grains 
as given by Greenberg and Hong (1973), are in the elongated core-mantle 
particles with silicate (or metallic oxide) cores of typically 0,05 pm radius 
surrounded by mantles of mean thickness about O.t pm composed of complex 
molecular mixtures of predominatiy oxygen, carbon, and nitrogen with hydrogen 
(Greenberg, 1978). The mantles may achieve a mean thickness as large as about 
0.15 pm (total grain diameter « 0.4 pm) at the time when coagulation of the 
duat leads to the formation of the comet (Greenuerg, 1979). Thus, the relative 



f 



volume of mamtle-to-core material is of the order of (^cl) 3 or about 60; 1, The 
index of refraction of the silicate core is approximately m «* 1.6 and the index 
of refraction of the mantle is about in » 1,3 as characteristic of complex icott i 
The absorption of the co fa and mantle material in the visible is quite low, 

‘The other constituents in. the bi-modn’l distribution are in a large number of 
vary small particles of not completely determined chemical composition, which 
however contribute at most a small amount both to the bulk (< IX) and the optical 
characteristics of the aggregated material (Hong and Greenberg, 1980), even if 
they axe predominantly made of graphite. This can also he seen by noting that the 
total absorption by such particles in the visual, in the interstellar medium, 
is small compared with the extinction by the core-mantle particles. 

Subsequent to breaking away from the comet, the debris will be subjected 
to both solar radiation and solar wind particles, This is o rather complex 
problem to treat in its entirety because one should follow the particles during* 
their entire lifetime 10 6 years) while their orbits bring them within varying 
distances from the sun, It is illustrative that at 1 AU the sputtering or 
erosion of such components of the aggregate as HO ice could occur at a rate 
of about 5A per year as a result of the solar wind (Umzerotti et al,, 1978), 

The leas volatile components like the silicates and the very largo molecules 
(which are not melted or evaporated at temperatures as high as 600 K) in the 
photo-processed grain mantles (Hagen et al,, 1979) would undoubtedly erode at 
a much slower rate. Thus comet debris particles of initially 100 pm radius 
could have parts of their outer material eroded away in varying amounts tc 
depths up to ^ SO pm leaving an open outer structure containing mostly the 
silicate cores and the nonvolatile mantle molecules. This tangle of elongated 
silicate particles as exposed in the outer portion of the comet debris particle 
or as imbedded in various degrees within the grain, mantle material, we envision 
as a bird's-nest type of structure. We note here that not only erosive effects 
can be produced by the solar wind particles but also as a processing of the 
debris matter analogous to the photo-processing of interstellar grains. The low 
density of the agglomerates as well bb of the fluffy particles studied by 
Giese et al. (1978) is apparently not consistent with the densities inferred 
by Le Sergeant and Lamy (1978) from lunar rock cratering but is consistent 
with the densities o- 1 g cm” 3 inferred by Nagel et al, (1975) and with the 
microroeteoroid detections by the Helios 1 spacecraft (GrCln et al,, 1980), 
where it was found that a large number of particles have densities below 
t g cm” 3 . In chapter 5 the optical properties of "bird's-nest" structures are 
compared to those inferred for the zodiacal- light particles. Clearly, the 
theory required to predict the exact scattering by such structures is far 
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Development of an "ice"- 
free outer crust 

3rin and Mendis (1979) 



No evolution 

Whipple (1978) 

Delsemme and Rud (1973) 


No change in grain chemistry 
during or after fragmentation 


Comet fragmentation 



in grain chemistry 
after fragmentation 


lErosion of 
volatiles 



Aggregates of loosly packed 
elongated silicates with 
inclusions of the smaller 

bare grains 

» 


Aggregates of elongated core- 
mantle grains with inclusions 
of the smaller bare grains 


Fig. 1.1 Successions of theories and assumptions leading to the 
concept of "bird's-nests". 
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beyond the scope of present techniques, Not only need one consider & large 
number of particles simultaneously, as in normal multiple-scattering theory, 
but one must also take into account that the embedded silicates are close 
enough to each other to require inclusion of near-field modifications. Faced 
’with this impossible situation we resort to the microwave analogue method 
described in chapter 2. Some of the targets discussed in chapter 4 are designed 
to clearly separate higher-order scattering effects from relatively easily 
predictable effects. The latter effects are included in the approximate method 
developed for investigation of light-scattering functions for some categories 
of inhomogeneous targets that fall outside the present laboratory constraints 
and especially for providing a first approximate understanding of this exotic 
field in light-scattering theory. The higher-order scattering effects are found 
to be negligible, to the comparisons with measurements in chapter 4 are quite 
reasonable. 
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Chapter 2 

THE MICRO-WAVE ANALOGUE METHOD 
AND 

THE LABORATORY 


2.1 Historical Background 

In I960, J. Mayo Greenberg constructed the first microwave analogue scattering 
facility to measure the extinction by single, irregular particles (see Lind, 
1966). Over the years, such i> inurements have contributed significantly to 
characterizing (size, chemical composition, core-mantle nature of) interstellar 
grains. The micro-wave laboratory under its director Dr. D.W. Schuerman is now 
a part of the Space Astronomy Laboratory (S.A.L.) which is a research-only 
organization of the State University of New York at Albany. Most of the work 
reported in the present study were made at this laboratory. 

The measurements for the three largest agglomerates were performed in the 
micro-wave facility of Bereich Extraterrestiche Physik, Ruhr Univeraitllt Bochum, 
German Federal Republic. In this laboratory scattering from cubes and other 
non-spherical particles, including the fluffy model for zodiacal-light particles, 
have been extensively investigated by R. Giese, R. Zerull and their co-workers. 

2.2 The Analogue Method 

In the general theory of scattering based on Maxwell's equations, all parameters 
involving the size of the particle are measured in units of the wavelength of 
the incident radiation. Thus, it is possible to scale the entire problem of 
visible, ultra-violet and infra-red scattering off micron and aub-micron 
particles to larger wavelengths and particle sizes. At S.A,L. the scale factor 
is about 10 s . Thus, apple-size particles are irradiated by micro-waves of 
X “3.18 cm. A particle of arbitrary shape is simulated by machining or molding 
a micro-wave target which is ~ 10 s times larger than the real particle. The 
target material is selected by requiring it to have, at the frequency of 
9.42 GHz, the same properties (dielectric constant, conductivity and size to 
wavelength ratio) as the actual particle has at visible, ultra-violet or infra- 
red wavelengths. 


v 
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2,3 Measured Quantities 

The scattering process is described by a 4 x 4 scattering (or F) matrix which 
operates on the intensity of the radiation as represented by the Stokes vector; 
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where k * 2r/X and R a is the distance between the scatterer and the receiving 
antenna. 


The matrix contains only six independent functions because it is assumed 
that the elements are averaged over a random distribution of similar particles. 
Since the input radiation, defined by the Stokes vector (1°, 1°, U ° , V°) , is, 
usually natural light, we set (1°, 1°, U° , V°) * (1°, 1°, 0; 0) with 

Q Q Q 1 i 1 2 

1^ » ■ I . From the form of the F matrix in equation (2.1), it is obvious 

that the Stokes vector of the scattered radiation must be of the form 

(J , J , 0, 0), so that, equation (2.1) can be written more economically as 
1 z 
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The measured quantities are I (0) , J (0) , 

11 12 

total ’-brightness 


•i 

and I (0). Given these values the 
22 


J (0) - I (8) + 21 (0) + I (0) 

tOt 11 12 22 


and degree ,of polarization 


I (0) - J (0) 

F(0) , - 

X (0) + 21 (0) + I (0) 
* 11 12 22 


(2.3) 


(2.4) 


can be determined. 
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2 .’4 Definition of Target Orientation 

The orientation of a particle is defined with respect to direction of the 
incident plane wave and of the measured scattered ra; /tion Figure 2.1 
represents an arbitrarily orientated cylinder whose symmetry axis makes on 
angle t|/ with the perpendicular X to the scattering plane k Q ** k g '(X"~Z) whose 
projection in this plane makes on angle x with the incident ^-vector. The 
angle between 7e and is the scattering angle 6. 



Fig. 2.1 Angular coordinates specifying the orientation of the symmetry 

axis op a cylindrical scatterer with respect to the direction 

of propagation Ji of the incident and Jc of the scattered 
o a 

radiation. 


In the laboratory a target is mounted in such a way that it may be 
rotated around two mutually perpendicular axes, the laboratory fixed 
azimuth axis X, and the target fixed tilt axis which sweeps in the scattering 
plane. For particles that do not possess an axis of symmetry, an arbitrarr 
axis is used as reference. 

2.5 The Microwave Laboratory • 

2.5.1 file mic rowave labor at ory at S.A.L . 

The target, in operation, is suspended in the center of the micro -wave- room 
by nylon monofilament lines, which are sufficiently thin to cause negligible 
scatter. The target is irradiated by micro-waves generated by a single screened 
dipole at the focus of a 51 cm diameter parabolic reflector located 9.84 m 
from the target. The dipole, which can be arbitrarily oriented in the focal 
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plane by means of a rotary joint, is fed through a waveguide by a 9.42 GHz 
stabilized oscillator. The resulting scattered radiation is measured by means 
of a square-law crystal detector fed by a similar 41 cm diameter antenna at 
5.18 m from the target. The signal ia chopped by a 1000 Hz square-law raodulator- 
■ -fed diode switch located between the detector and a hybrid junction. In the 
junction the scattered signal is mixed with a background compensation signal, 
made equal in amplitude to Che signal from the receiving antenna in the absence 
of tho target by means of a variable attenuator and made destructive by use of 
a variable phase shifter. The compensation wave is tapped fron,' the waveguide 
supplying the transmitting dipole by a 10 dB directional coupler. The length 
of the waveguide is choosen to minimize the optical path differences between 
background and compensation signals at small scattering angles, where direct 
radiation may be intercepted. A rotary joint on the waveguide permits the 
entire detection system to he mounted on a moving carriage and located at any 
scattering angle 0 from 0 to 170°. A second antenna and detection system is 
located directly in line with the transmitter*' target; it measures the amplitude 
and phase of the forward-scattered (0 “ 0) radiation. This antenna is located 
in a radiation Crap and is covered by a micro-wave absorber during the side- 
-scattering isaasurements. Further technical details may be found in Greenberg 
et al. (1967) and Wang and Greenberg (1978). 

As part of a continuous effort, increased experimental accuracy and speed 
has been achieved by digitalizing the data handling and target orientation 
mechanism in two phases during the measurement period. We name the versions 
used in this study A, B and C. 

2.5.2 Versio n. A 

The diagram in figure 2.2 represents the apparatus of version A. The output 
potential over the detector, which is proportional to the intensity of the 
incident radiation, is amplified up to 10 7 tiroes by a lock-in amplifier unit 
linearly feeding the ordinate on an X-Y plotter. The abscissa is connected to 
a potentiometer which senses the azimuthal angle (x) driven by an analogue 
motor. The motor is controled by a manually operated potentiometer through a 
feed-back mechanism using the X"detector. A similar system controls the 
elevation and the tilt angle. For a detailed description of the orientation 
mechanism we refer to Lind (1966). 



Fig. 2.2 Hicrovave circuitry version A, antenna diameters and ranges for angular 
distribution studies, at Space Astronomy Laboratory, State University of 
New York at Albany, Albany, New York. 
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2.5.3 Ogetat i^ - grocedure8 x- vers ion^A 

The detector feeding a lock-in amplifier, described by Brophy (1977) or in 
any basic textbook in electronic engineering, is intermittently screened out. 

In the process of amplifying, the noise signal thus detected is subtracted. 

•The square wave feeding the diode switch is used ns reference signal. First, 
the phase is adjusted so that the output signal vanishes; we are then 90° 
from the correct phase if the reference wave is symmetric. At 180 phase shift 
the true signal 4 s subtracted from the noise; the symmetry of the square wave 
is adjusted so that the output U(180°) ** -U(0). At the beginning and end of a 
set of measurements the system is calibrated by detecting the known scatter of 
a standard target. Around the selected scattering angle 0 ** 45°, the phase 
function for the lucite sphere with circumference to wavelength ratio 4.978 
is relatively feature less and the background is low. In the absence of the 
target, the residual signal is brought to less than 0.1 pV using the variable 
phase shifter and attenuator on the compensation waveguide. The target is 
placed in the beam and the height controlled through a telescope. One of 13 
linear scales amplifying from 10 to IQ 7 times were selected. This one scale was 
used as the target was rotated n full 180° around the azimuth & is. As the 
target rotates small oscillations in the detected signal, caused by inhorao- 
geneties in the target material, average out. Subsequently, the scatter from 
targets with unknown optical properties can be measured. The background 
compensation procedure is repeat d at eacn new scattering angle, and the 
apparatus is re-calibrated fo* each combination of polarizations. Intermediate 
calibrations were performed during long sets of measurements, especially after 
changes in the external conditions such as sunrise, sunset or thunderstorms. 

The data reduction from the graphs were performed manually (~ 10*' data points 
pen target). The data presented for targets U41R and U41a was accumulated in 
this way. The measurements for targets TMR, T21R, ThIR, T4lA and C41 using 
this apparatus were in general agreement with the corresponding data using 
versions B or C of the laboratory. The later data is used in this study. 

2.5.3 App aratu s and operating pggggdureg ± _vgrsion_B 

The analogue motors in the orientation mechanism were replaced by step-motors 
orienting the target with a resolution of 0.1°. The tilt mechanism, originally 
used for very elongated targets, was redesigned and the range increased from 
< 180° to 180° by use of a constant tension spring-motor counteracting the pull 
of the step-motor. The range of the azimuth drive was increased to 400°. Two 
sets of micro-wave absorbers were added to screen out the orientation mechanism. 
A PDP-11/03 DICLAB minicomputer was installed to control the step-motors, 
reading of the Xock-in amplifier, and for data reduction. A laser for re- 
alignment of the receiving antenna at each scattering angle was mounted on the 
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the dashed lines are the nylon e*onn filament string# supporting and 
orienting the target* 

figure codes 

a - spherical target 
b ** target orientation mechanism 
c - to the hoist 

d «• transmitting entenrvi 
c - receiving antenna 
t *> background compensation waveguide 
g * micro- wave absorbers 
h * wini’-cowputcr 

x » power supply and computer interface with the step motors 
j - detector electronics on * while rack 
k - telescope for target alignment 
l - laser for antenna alignment 
ra - target rack 
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mobile carriage. The computer was programmed to step the target through any 
preselected set of orientation. At each orientation 21 measurements were 
collected, and if the mean fell outside 15-85% of the full-scale value at the 
selected gain the measurement was rejected and new measurements were taken at 
this orientation when the appropriate scale (requested by the computer) had 
been set. If not rejected, data was sampled at each of these orientations 
until the standard deviation of their mean was less thon any preselected value 
or a maximum of 10'* measurements was reached. The average and the number of 
meaourements were stored on a floppy disc along with a flag indicating if the 
preselected standard deviation of the meat as reached. Using this apparatus 
targets T11R, T21R, T4IR, T41A and IC la required an Average of 250-300 hours 
each for a full set of measurements, 

2.5.5 Apparatus _and p< opera .ting procedu re s, ..version C 

A second lock-in amplifier controlled by the computer and with maximum gain 
of 10 e was installed. The manually controlled amplifier was still used for 
nulling. This and other refinements have gradually increased the rate of 
measurements to ~ 45 hour/particle, 

2.5.6 Errors 

The uncancelled background signal X ^ < 0.02 interferes with the scattered signal 
I g , The resulting intensity X is given by 

vr - vt b +\/i h cos A 

where A is the phase difference. If the intensity I g incident on the detector 
causes a potential difference of 10 yV. The potential 8.1 < Uj. £ 12.2 in yV 
or in absolute units 0,162 £ i 0.244 results. How'ver, the phase of the 
background signal may vary rapidly (if not rapidly variable it can be cancelled). 
If A is truly random over the set of measurements, .at each orientation the 
intensities add directly; thus, Uj. ■ 10.1 yV, In the case of the "bird’s-nests" 
at higher scattering angles, the phase of the scattered wave takes all values 
as the target is stepped through the 272 orientations, therefore the maximum 
error of the mean is given by (\4^ + v/z^) 2 - I g at 0 * 0 and approaches 
= 0.02 at higher scattering angles. Drift in the gain may reach 10%. For 
this reason, the system was re-calibrated against the standard sphere inter- 
mittently during the course of measurements. The (generally) small cross- 
-polarization comp' nents vanish for spheres, the average of the calibration 
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constants for the 11 and 22 components were substituted. The agreement between 
theory and measurements for scattering by spheroidB using verions B and C of 
the apparatus was found to be very good (Schuerman eU al., in preparation). 

This was still true when the cross-polarization terms were calibrated a6 in 
tnis study. 

2,5.7 Thej^ijcj^j^^ji_laborntory_nt_^^ 

The laboratory operates at 35 GHz and is intended for side-scattering only. 

Both transmitting and receiving antenna are of the horn type. The mobile 
receiving antenna, mounted on a track, stops automatically every 5° as it 
moves from 0 m 155° to 50°, In this interval, there is no compensation for 
the background. At other scattering angles the target is removed from the beam 
and the background cancelled in a manner similar to that at S.A.L.. The 
scattered signal is integrated electronically as the target is rotated 360° 
in the scattering plane. The gain is calibrated through the known scatter of 
a standad sphere at one given scattering angle. For a more complete description, 
see Zerull (1973). 

2.6 Target Manufacturing 

2.6.1 ’’Bi rd 1 3-nests" 

In space, the "bird’s-nest" type particle problem involves three media, vacuum 
(jn » 1) , "ice" ( m & 1.3) and silicates ( m « 1.6). This is modelled in the 
laboratory-, using the four media; air (/» w 1), eccofoam (m « 1), eccofoam PS 1.7 
(m & 1.3) and lucite (m & 1,6). In all cases the imaginary part is small. The 
eccofoam is used to fix the relative positions of the individual scatterers and 
is similar to eccofoam PS 1,02, which is optically almost identical with air 
except at small scattering angles. The match of refractive indices is within 
the uncertainty for the cosmic materials. The "bird's-nests" were built in 
modules (see Table 1) which can be combined arbitrarily, Lucite, also known as 
plexiglas, is an acrylic that is commercially available in long rods. The 
cylinders were cut to size from appropriate diameter rods. The indices of 
refraction ( m * 1.618 - 0.005-£) of two samples cut out of an 1.5 inches diameter 
cylinder, both made to fit snugly into a waveguide only differing in length, 
were derived from two sets of standing wave patterns arising inside the wave- 
guide containing a respective sample using the slotted waveguide method of 
Roberts and van Hippel (1954) and further developed by Wang (1968). The eccofoam 
PS is a low loss polystyrene manufactured by Emerson & Cuming Inc. The dielectric 
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constant e, being closely related to the density of the target, is adjusted 
by pre-heating, and thus controlled expansion of the target material, Eccofoam 
PS 1.02 has e « 1.02. In eccofoam PS 1.7, microspheres (36 DS) are mixed to 
.the polystyrene to achieve e - 1.7 to an accuracy of ± 0.02 (Emerson & 

Cuming Inc., pers. com.). The dissipation factors tan 6 are smaller than 
0.0002 and 0.0001 respectively. The resulting refractive indexes, 


m 


\f s 


e \/e 2 *■ (e tan 6) 2 


- t 


V 4 


+ \/e 2 -T (e tan 6) 2 * 
2 


are « 1,01 and 1.304 ± 0.003 respectively. The eccofoam PS 1.7 was sawed very 
slowly (less than 1 inch/min. ) and cooled with running water to avoid further 
expansion due to friction heat, The resulting square rods (1 x 1 x 5.5 inches) 
were machined circular in a lathe and cut to 1.0 or 0.8 inch diameter cylinders. 
Each cylinder was cut in two, and a concentric hole drilled (se« Fig. 2 , 7 ) 
to accomodate the lucite cylinder. The two eccofoam parts were glued together 
using a paper glue with index of refraction close to 1.3. The resulting 
concentric cylinders, with length-to-diameter ratio of both core and muntle 
2:1, were exposed to a 0,4 mm Hg vacuum for more than an hour to evaporate any 
remaining water and to dry out the glue. 



Fig. 2.7 Model of an interstellar core-mantle grain. Each eccofoam ( M icy n 
material at visual wavelengths) cylinder is cut in two and a 
concentric and coaxial hole drilled Co accomodate the lucite 
(modelling silicate type material) cylinder. The two eccofoam 
parts were glued together using a paper glue with index of 
refraction close to that of the eccofoam (n«1,3), 
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Fig. 2.8 One of the spherical eccofoam shells pictured as part of the 
target positioning mechanism. The nylon monofilament string 
connecting to the tilt motor windes on a machined circular track 
and reinforces the bound between the Upper and lower eccofoam 
hemispheres by first passing diametrically across the sphere then 
following the machined track 360° and passing back diametrically 
across the sphere again. The nylon (fish-) line and the eccofoam 
are almost invisible at the microwave part of the electro-magnetic 
spectrum. 
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The hollow eccofoam spheres (m « I define the shape and size of the 
agglomeration and ere part of the suspension mechanism. Homogeneous spheres of 
this material are commercially available. Using an arrangement similar to that 
of a lathe, a well defined track was carved, out for the tilt string. Then the 
, spheres were cut approximately along a great circle and the halves carved out 
with a razor blade to avoid crushing the materiel. The spherical boundary was 
purposely made rough to prevent the cylinders sliding. The halves were glued 
back together and an opening cut out. The volume was measured by inserting a 
plastic bag and filling it with water to obtain the volume (or weight). The 
removed material was saved to be inserted between cylinders. The hole passing 
the string connecting to the hoist was ruinforced by small pieces of cardboard. 
The string connecting the tilt mechanism also reinforces the bond between the 
hemispheres as pictured in Fig. 2,8. The spherical cavity wan then filled with 
the desired mixture of cylii'.ders and filling material (eccofoam) such that the 
number density of cylinders was constant in layers. 

Some of the smaller "bird’s-ney'ts" containing only lucite cylinders were 
made by inserting the cylinders into an eccofoam sphere thereby keeping the 
number density constant in concentric spherical layers Two, orthogonal annular 
rings were fit around the sphere (no glue) . On one the track for the tilt 
mechanism string was machined; Che other was used to attach the strings to 
the hoist. The two ,! bird ! s-nests” containing aligned cylinders were also made 
from solid eccofoam spheres but these were cut in hemspheres, and the cylinders 
were inserted perpendicularly to the plane surfaces to facilitate alignment. 

As these targets were never tilted the strings attached to the hoist were 
glued directly to the ec»„o’oam. 

2.6.2 Sin g l e cyl i nde r s 

The single circular homogeneous or concentric cylinders, as described in the 
previous section, were suspended by two strings glued directly to the cylinder 
on diametrically opposite points equidistant from the two ends. The two strings 
used to control the target tilt-angle were glued on opposing cylinder edges, 
such that the projection of the total of four points of attachments on a plane 
perpendicular to the symmetry axis are corners of a square. 

2.7 Simulation of an Optically Thin Cloud of Randomly Orientated Scatterers 

An optically thin cloud of randomly orientated particles is simulated by 
averaging the signals received from a few hundred uniformly distributed particle 
orientations. This set of orientations for a particle having axial symmetry 


is obtained by pointing the symmetry axis of the particle to 272 positions almost 
uniformly distributed over an imaginary sphere centered on the particle. Perfect 
uniformity cannot be achieved - a fact known to the ancient Greeks who discovered 
that only 4, 6, 8, 12 or 20 points could be so distributed over a sphere. 

However, by distributing points over a pentagon in an almost uniform manner, 
one can combine 12 such pentagons to form a dodecahedron - one of the five regular 
figures which can be inscribed in a sphere. In this sense, we obtain the closest 
approximation possible to uniformity. 

2.7.1 ygifQ£m_distribution_o£_alJL_orientations_of. the single cylind er s 

The single cylinders C41 and 12 have both an axios of rotational symmetry and 
a plane of symmetry perpendicular to this axis. The cylinders were orientated 
so that their axis of symmetry was successively pointing to 1/8 of the 272 
orientations determined as described above. The 34 orientations were chosen 
as described in appendix A for axial symmetric particles with a plane of 
symmetry perpendicular to the symmetry axis. Also 8 fixed orientations including 
alignment of the symmetry axis with the incident &, 2? and $ -vectors were 
investigated; but the results are excluded from the present study. 

2.7.2 Unifo rm distrib u tion^ ..orxentacigns o£ a. ''bi rd 1 s-nest” 

The "bird's-nest” type dust grains are asymmetrical. Thus only the symmetry 
relation a defined in appendix A prevails. All possible orientations using 
the laboratory orientation mechanism are independent as long as the tilt axis 
does not coincide with the bisectrix of the complement to 6. We will first 
consider the orientation of "bird ' s-nests" consisting of cylinders aligned for 
perpendicular incidence. To keep the cylinder axes perpendicular to the incident 
radiation we- are only free to rotate the target around the azimuth axis. Let 
X “ 0 be an arbitrary orientaciou around the azimuth axis. The target is then 
orientated to' 360°- Ax with increments of Ax. At each orientation, measurements 
are performed until the standard deviation of their mean is less than a pre- 
-selected value (with a minimum of 21 and maximum of 10,000 measurements). As 
can be seen in Fig. 4.2 the scattered intensity oscillates heavily as the 
target is rotated, as a result of interference of the contributions from 
individual scatterers. The periodicity p of the oscillations is dependent on 
the scattering angle. This phenomenon is discussed in appendix B. To avoid 
systematic errors we choose Ax « p. This condition leads to several thousand 
orientations at the higher scattering angles. If the "bird's-nest" consists 
of randomly orientated cylinders we may also rotate around the tilt axis, or 
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•to 10* - 10 7 orientations. At an average rate of 10 sec per orientation, the 
tine required is of the order of a year or more. Instead the 272 orientations 
resulting from a ctaMnation of 12 pentagons were used. 

2.8 The Targets Selected 

Of the 17 targets used in the measurements, eleven represent agglomerates 
of, or single, homogeneous silicate cylinders. The remaining six targets 
represent agglomerates of, or single, silicate-core - ice-mantle cylinders. 

It is convenient to refer to 21 separate sets of measurements as represented 
in the fold out page, because sorn targets have been investigated at two wave- 
lengths or at differing orientations. In the group of ten targets nine sets of 
measurements (BAIR, U41R, C41R, B21R, B11R, U41R, C41A, P1 1 and T41A, where the 
two digits stands for the elongation of individual cylinders and the final 
letters A for aligned and R for random) and one (12R in the group of six targets) 
were selected to illuminate the multiple scattering problem and provide tests 
for the theory. Target 12 represents a single interstellar grain, and its 
scattering characteristics are therefore of interest for other astronomical 
applications as well. The rest of these scatterers do not represent likely 
cosmic grains at visual wavelengths, with the possible exception of T41A. 

All the other pure silicate "bird's-nests", model cometary debris from which 
the volatiles have erroded before or after fragmentation. The silicate-core- 
-ice-mantle group of targets represent pieces from comets that have remained 
premordial . 

The agglomerates correspond to the 1.5 - 6 pm diameter range at X * 5000 A 
and our biggest model for actual cometary debris corresponds to a 3 pm particle. 
This is dictated by the laboratory constraints. Actual comet debris, however, 
have a rather wide size distribution and the optically most important zodiacal 
light particles are thought to be in the 10 - 100 pm size range (Weinberg and 
Sparrow, 1978). For future interpretation of scattering from larger dust 
par, tides an approximate theory is developed in chapter 3. 


Chapter 3 


SCATTERING APPROXIMATIONS FOR 
SOME INHOMOGENEOUS PARTICLES 
OR SCATTERING SYSTEMS 

3 # f Scattering at Lour Wavelengths Compared to the Inhomopeneties of Thei;. 
Mutual Distances 

We may divide the inhomogeneous target material into volume elements that are 
small compared to the wavelength X of the incident radiation. If the material 
i6 randomly inhomogeneous there is a wavelength such that the integrated 
polarizability over the volume elements is uniform. At this and larger wave- 
lengths the particle acts like a homogeneous scatterer with the index of 
refraction m given by 

m « 1 + 2v P , (3.1) 

where P is the induced dipole moment per unit volume. 

O fl mi w— i. * a 

iuc largct natcnai wnen m re i 

When X is made short enough some volume elements contain individual molecules 
while other are vacuum. If the mutual distances between volume elements with 
m j* 1 are > X each volume acts like an individual scatterer. The extinction 
and dispersion in a dilute medium of N randomly spaced scatterers with forward 
scattering amplitude 5(0) is discussed by van de Hulst (1957, p. 32). It is 
shown that a wave travelling through the medium suffers a phase lag and 
attenuation, that can be represented by the complex index of refraction m of 
a homogeneous material, 

• m « 1 - i 5(0) — . (3.2) 

k 3 

Let z - 0 be the boundary of the target material. Then we may write the 
incident electric field 

E Q - e Ut (3.3) 

and the field incident from the negative z direction on Che particle with 
coordinate z/ ., 2 . in the scattering plane as 
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Combining eq. 3.3 and 3.4, 

E . - B . 

0 % 0 

The field acactered from the target to the point y,z in the far field is 
given by 


( y ~ y i\ 

E q “^arctan— J 


y~Vj\ m + iwt 


ik\{y-yj) z + (z-z.) 2 


or in units of the radiation incident on the acatcerer, 
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Insertation of eq. 3.5 into 3.7 gives 
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The total scattered intensity from N particles is simply 


Ky,z) - ^ . 


In the forward direction y**u • thus eq. 3.8 become: 

% 


-ikm 


E • s( 0) — E 

8 iUz-z.) 0 


(3.4) 
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(3.7) 


(3.8) 


(3.9) 


(3.10) 


so that the amplitude but not the phase of the forward scattered wave depends 
on the particle location. Thus the scattered radiation is strongly peaked in 
the forward direction because of interference. Note that in the limit where 



Che material of infinite dimensions Is homogeneous 1(0), 0 f 0 vanishes. If 
the Rayleigh approximation S(0) • ik 3 a lor the forward scattering by the volume 
elements is substituted the mean index of refraction can be written from eg, 3.2 


m * 1 *► 2ifotf (3,11) 

where a is the polarizability of the particle. The assumption made in deriving 
eq. 3.2 is that the system of scatterers is dilute, which impliea that m * 1, 

In this limit the Loren tz-Lorenz formula 

. 4tto# » 2$LzJ l (3.12) 

m 2 + 2 

can be rewritten as 

4xcdV w m 2 - 1 w 2 (m - 1) (3.13) 

so that the two formulae are equivalent for m « 1. This relationship is pointed 
out by van de Hulst (1957). 

So far we have discusived the extreme cases where X > the dimensions of the 
individual scatterers and further more their mutual distances are either X 
or > X. The discussion and the formulae in this section apply to any dilute 
system of scatterers separated by distances large compared to X and to the 
particle dimensions with the appropriate solution for £>(0). 

3.3 Particles Very Large Compared to the Wavelength 

The condition in the previous section was that the distance between scatterers 
» X. Wits implies that the formula derived in that section only apply to 
scattering systems » X. In this particular case the incident plane wavefront 
of infinite extent can be thought of as consisting of separate rays that persue 
their own path. The condition is that the ray of length £ must have a width of 
the order of VlX corresponding to the first few Fresnel zones. The rays passing 
along the particle form a plane wave from which a part in the form and size of 
the geometrical cross section of the particle is missing. This incomplete wave- 
front by Huygens' principle produces the Fraunhofer diffraction pattern at 
small scattering angles (see van de Hulst 1957* p. 24, or any textbook in basic 
scattering theory). The Fraunhofer diffraction is unpolarized, and by definition 
independent of the internal structure of the particles. Furthermore, the 
restrictions on target shape are not very rigorous. At larger scattering angles 
the scattering resulting from eq. 3.9 where the summation is made over the whole 
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ensemble, dominates. An interesting result is known as Babinet's principle 
(van de Hulst, 1957). We will make the classical set of imaginary experiments; 

(1) we replace the particle by a disc that totally blocks an area of the wave- 
front of the same shape and size as the geometrical cross section. The diffrac- 
tion observed is identical with that of the particle blocking the wavefront. 

(2) If instead the whole wavefront is covered, except for a hole of the same 
geometry bb the particle cross section, n similar diffraction pattern is seen, 
but the electric fields are of opposite signs. (3) Obviously, if we were to 
combine the two diffraction patterns (1) and (2), they cancel each other by 
destructive interference, In (3) we do not block any part of the wave one 
therefore the wavefront must remain undisturbed. Let a be a typical dimension 
for the particle (for spheres we let a be the radius). Then similarly, if we 
place a largi particle with the index of refraction m < 1 + in the homo- 
geneous wavefront, the diffraction pattern (1) is partially cancelled. The 
restriction on m means that the phase difference between light rays that havv 
passed aside the particle and the transmitted wave is small. 


The resulting scattering pattern for the large scatterers discussed in 

this section is strongly peaked in the forward direction where the Fraunhofer 

diffraction pattern (1) dominates. If the condition m 1 + is fulfilled, 

the diffraction pattern has the same shape as in (1) but is lower in magnitude. 

In other directions the scattering given by eq. 3.9 dominates. The side- 

-scattering depends strongly on particle shape and size and goes generally 

through a large number of oscillations as a function of 6. As an example, an 

ensemble of scatterers in the shape of a hexahedron aligned along the axes 

X,X,Z, defined in section 2.4 will scatter in phase at 6 <* 0, 1(0) - (,£ E.( 0)) 2 . 

If the target contains N truly randomly located scatterers and the thickness is 

n kj 1 , where n is a lar^e integral number and \ ff/ is the wavelength inside the 

target, lim(r(l80°)) * E (180°) « .1 I. (180°), It the dimensions along the 
tf-Ho l N o 

X- and Z-axes each are mupliples of A . then lim (1(90 )) ■ £ I. (90 ). If m 

m im 1 1* 

is 'isotropic (randomly oriented scatterers) eq. 3.9 is indpenendent of the state 
of polarization of the scattered radiation. Therefore the degree of polarization 
is smell where the Franhofer diffraction domintates and approaches that for 
individual scatterers as 0 increases. 


3.4 Expeeted Effects on the Scattering Pattern as the Target Size Decreases 

The Fraunhofer diffraction pattern will, extend to higher scattering angles 
and decrease in relative Local intensity when the target size decreases. The 
interference pattern (eq, 3.9) remo f ne constructive father out and the number 
of oscillations in the scattered intensity as a function of 0 decreases, As 
the number of individual scatfererc decreases, the angular distribution 
becomes increasingly dependent on the relative location of scattcrers, because 
the phase differences between contributions from individual ocatccrors 
fluctuates statistically around the asyi„j ^ntic mean value at as different 
but always random distributions of the finite sample of acatferers are chosen. 
This effect is observed in the laboratory as oscillations in the scattered 
signal when the ensemble is rotated at a fix scattering angle. This effect 
and it dependence on 0 is discussed in appendix B. 

3.5 Expected Effects on the Scattering Pattern as the Distance between 
Scatterers Decreases 

The components of the polarizability vector joining individual scatterers, and 
multiply scattered radiation, will increase in importance when the distances 
between scatcerers decrease. These two effects car. not easel y be evaluated 
partially because th® particles are not in each other tarfielde any longer. 

We resort to the micro-wave laboratory, 

3.6 The Single-Scattering Approximation 

Both single- and secondary-scattering computations were performed for T41A 
consisting of 125 randomly located 4; 1 -elongation silicate cylinders all 
aligned perpendicular to the scattering plane, each with the circuraference- 
-to-wavelength ratio or size parameter x * 0,47 in a spherical ensemble 
with size parameter X n 10.04. For each cylinder the scattered intensity 
and phase respectively were computed separately. This was dictated by the 
limitation of core size available in the computer. 

3.6.1 Intensity 

For the finite cylinders with x *= 0.47 the particles are too small to provide 
reliable experimental data as individuals. Also the infinite, cylinder 
computations are not readily normalized (see appendix D) because the intensity 
at a distance r in the far zone is proportional to r"" 1 rather than X** z as for 
finite targets. Therefore, we have substituted appropriately normalized 
spheroid calculations. Each cylinder is replaced by an equal volume spheroid 
whose scattering is computed using the Rayleigh approximation. 


3.6.2 Phage 

The phase la defined with respect to on imaginary wave deflected at the target 
site in vacuum. The optical path of radition scattered off a cylinder A 03, a, 3) 
in the "bird's-nest" of radius R and with moon index of refraction m ^ * m' « irn 

is 

0 * UijA + i7i n ) ( 1 ~ n’) + iS^A + 4'tf,, (3.14) 

sec Pig. 3.1 . 

Let r *» \//? ? « i^sin 2 ((i) 

a * Pcos(0) cos(r) 

y ® Z)cos(p) sin(a) (3.15) 

Then} 

e^A ** Vr 2 - j / 2 + x 

_ ■ ■ « > - 

Ae^ * vr 2 - (//cos (0) + «sin(0)) 2 - «cos(0) - ysin(0) 

t)"^A * Vh i ** y 2 + a 

4^ * '/t? 2 - (j/cos(0) + .rsinCfi)) 2 ~ «coa(0) - ysin(O) (3.16) 

The resulting phase, where ph^ s,n the phase shift over the cylinder, is 

PH » CO - W)~ + p/i 4 (3.17) 

Note thatOT 1 in eq. 3.14 formally results from secondary scattering. 

3.6.3 Some j)rediction£_by_the_single_ncattering_aggroximatian 

The maximum phase shift in the 0*0 direction over the ensemble with radius 
R is given by 

Pfl ma X - 2 '%- ,)X “ ,6 ¥ P (3 - 16) 

where the dipole moment P per unit volume is proportional to NR~ 3 and N is the 

number of identical agglomerated particles, thus PR «NR~ Z and the inter- 
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Fig. 3.1 The single scattering geometry, A is the location of the 
scattering cylinder with spherical coordinates D,b t,f3. 

Aw ^ is the perpendicular distance between A and the incident 
wavefront tangent to the "bir-d-nest" with radius /?. Aw q is 
the perpendicular distance between A and an imaginary scattered 
wavefront tanget to the "bird' s-nest". 
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ference pattern becomes increasingly constructive as the ensemble expands or 
N decreases. In the limit , eq. 3.14 becomes 0 ^ » W^A + AW Q , The position 

(D,a,3) of a cylinder can be written (xdf,u,|3). Then 

WjA " ft + Xdcos(R)eos(«) 

Aw q ■ R ~ Xdcoa((J)cos(ct -6), (3.19) 

and 

PH » Xdcos(p) (cos (a) - cos (ot - 0))~ + ♦ (3,20) 

It is seen that in this approximation the dependence of PH on a and (3 vanishes 
at 0 » 0, where the interference pattern depends only on ph , the phase shift 
over single cylinders. Eq. 3.20 can be rewritten 

PH - ph £ - - ~~ cos(3)sin(ct - |*)sin(-|), (3.21) 

where it is seen that PR ^ and therefore the phase differences and their 
dependence on a and 0 are proportional to X and increase with 9. 


3. 7 Secondary Scattering Approximation 


3.7.1 Intens it y 

In the previous section each cylinder was assumed to scatter the light incident 

on the agglomerate independently of all other cylinders. However, some of the 

light may be absorbed or scattered by other cylinders. This is the phenomenon 

of extinction represented by the imaginary part of the index or refraction. 

The thus scattered light will irradiate the cylinders from directions deviating 

from that of light incident on the agglomerate as a whole. However, we assume 

that the cylinders '’'look" infinite at a distance comparable to their length so 

chat the cylinders are in effect irradiated at perpendicular incidence, Let 
% 

I ^ be the incident intensity on cylinder /l(i)^,a^,(3^) , then the intensity J ^ 
scattered of A and incident on cylinder B(D^,a^,^^) approximated by the far- 
-field solution is the cylindrical wave 


■Ir 


I A P(4>)(£ - l^sin(^) - fl B sin(# B )|) 


( 3 - 22 ) 


k \ jD^cos 2 (g^) + cosMgg) - 2 ^ 006 ( 1 ^) D b cos(0 B )cos(a^ - a g ) 


where l is the length of the cylinders and ip is the angle between the direction 
of the incident radiation and the direction to B. The parenthesis in the 
numerator denotes the length of B illuminated by A (the perpendicular projection 
of ,4 on B) if this is 0 or negative no secondary scattering is assumed to 
occur. The J element in JF^q)) varies only by 13% between tp = 0° and 180° (Fig. 
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PROJECTION IN THE 
XZ -PLANE 



PROJECTION IN THE 
VZ -PLANE 


Tig. 3.2 The geometry for the secondary scattering approximation. The 
light is first scattered by A and then by B. In the secondary 
scattering approximation is the length of B illuminated 
by A . 
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4.7) and to save computer core ^(p) is approximated by the mean value 0.095 
at all (p's. Also the phase function for cylinder B is approximated by this 
value. The intensity of the secondary scattering in 1^ is thus isotropic at 
the far field, at which the dimension of the "bird's-nest 1 ' can be neglected 
compared to the distance p to the observer. 'The intensity of the radiation 
first scattered off A and then off B is, 



— - {l- |^sin(^) - D g sin(3 5 ) I) 
k p » ~ 


(3.23) 


3. 7-2 Phase 

The. optical path is given by 

AB - m' \J ZTjeos 2 ^) + cos 2 (P s ) - 20^006(^^)003(^)008(0^ - a B ) 

(3.24) 


where m' ia the real part of the index of refraction as determined by eq. 
3.1 or 3.2. The total phase shift is then; 

+ Ad + Bw„ 

■ PH ■ Zv + ph^ + phg 


(3.25) 


where ph ^ =■ phg m ph is the phase shift due to the scattering by the cylinders. 
In the Rayieigh approximation ph * 0° , while for spheres of the same size- 
parameter as the cylinders, Mie theory predicts ph close to 1.7° in I at all 
0’s. ph ** 1.7° is used in these computations. 






Chapter 4 


EXPERIMENTAL RESULTS 


VERIFICATION OF THE ''INTERFERENCE" THEORY 


4 . 1 Agglomerates of Homogeneous Cylinders 
4 .'1 . 1 ^£_single_scattering_a£groxiinati.oii_a£2li£d_ t to_T4 I^A 

T4tA: 125 x - 0.47, 4:1 -elongation homogeneous cylinders aligned perpendicular 

to the scattering plane within an X “ 10.04 sphere, 

1 In the single scattering approximation the scattering from each cylinder is 
replaced by the Rayleigh approximation for an equal volume V prolate spheroid 
with the same elongation as the cylinders. The polarizability a ^ along the major 
axis j m 1 and the minor axis j “ 2 is given by 


m 2 -1 


(4*1) 


Let e 2 ■ 1 - ( b/a ) 2 where a and b are the length of the semi-major and -minor 
axes respectivly, then for the depolarization factors we have 


1 — - (-1 + J~ ln-i-±-£) 

e 2 2e 1 - e 


(4.2) 


L + 2L - 1 

X, 2 


(4.3) 


Using hfa « 1/4, it follows that: 


L - 0.0754, 

L * 0.4623. 
z 


(4.4) 


For details we refere to van de Hulst (1957, p. 70). Using the index of refraction 
(m ■» 1.618 - 0.005i) obtained from the larger casted cylinder (see sec. 2.6.1, p,19) 
in eq. 4.1 wa obtain: 
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a - 0.02596 - 0. 00023*. 

1 ' 

a « 0.01666 - 0.00010*. 

2 


(4.5) 


* The scattered fi'/ld entering in eq. 3.9 is obtained from the scattering tensor: 
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a cos(6) 0 

2 
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-ik 3 
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s 
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1 _ 
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(4.6) 




The cylinders are embedded in the expanded polystyrene matrix material with 

index of refraction m (assumed real). The mean index of refraction m of the 

o eq 

ensemble used to compute the optical path of light travelling inside the "bird’s- 
nest" is given by eq. 3,1 and may be written 


ro„„ m m + 2 irP. 
eq o 


(4.7) 


The dipole moments P pet unit volume for the electric vector parallel (11) and 
perpendicular (22) to the aligned particle axes are given by 


P - N 

n 


P ■ N 
22 




4ir 


V 


(l +-M 

\ 2 m 2 - 1 ) 


4tt 


V' 1 

ene 


Z ” 1 

ene 


(4.8) 


where V is the volume of the ensemble. 

esTZo 

Using the computed values of P for radiation polarized parallel and 
perpendicular to the spheroid axis, we finde 


m f 

* m 

•f 

0.060, 

m n » 

5 • 10- 

XI 

0 



n 

m* 

« m 


0.038, 

m n w 

2 • 10“' 

22 

0 



22 


(4.9) 


H8 expected, m n « m '- 1 so that the ensemble material acts like a dielectric. 
The scattered intensity given by eq. 3.9 is dependent on the orientation of the 
ensemble (defined by the angle x, see sec 2.4). The details of this scattering 


i 
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pattern depends on the exact location of the individual cylinders. In the 

theoretical target the cylinders were distributed using a random number 

generator, whereas in the experimental target each particle was inserted 

by hand in as close an approximation to randomness as possible. The resulting 

’X -dependence given by the theory vFig. 4,1) and the experiment (Fig, 4.2) 

are therefore seen to differ in detail hut to show the same over-all 

characteristics. The run of the standard deviation o . .(0) of the scattered 

tQ 

radiation J. .(0) , normflized to J, .(0) is a measure of the isotropy of the 

to %o 

optical properties of the scatterer. The ratio a . ,(0)/J. .(8) is independent 

to to 

of the gain of the apparatus and is likely to contain extensive information 
about the infernal structure and shape of particles. The predicted and 
measured run of with 0 are compared in Figs. 4.9 (,i m 0 m 1) and 4.10 

(t>j*2). To simulate a cloud of similar agglomerates, the scattered intensity 
at each scattering angle was averaged over 72 orientations. The resulting 
angular distribution from both experiment and theory is shown in Figs. 4.3 
and 4.4 for radiution polarized along the cylinder axes (perpendicular to the 
scattering plane) and perpendicular to the cylinder axes (in the scattering 
plane) respectivi} . The Fraunhofer diffraction pattern which is confined to 
the forward scattering hemisphere, and expected to dominate for larger 
agglomerates at small scattering angles is still seen and well represented by 
Mie-fheory, although the sice of the ensemble is as small as X ■ 10.04. 


If we replace the finite cylinders by equal length portions of infinite 
cylinders in the theoretical solution, the magnitude of the scattered intensity 
is increased by a factor w 3 (Figs, 4.5 and 4.6) because the scattering by an 
infinite cylinder it. concentrated precisely in the scattering cone while the 
scattering by a finite ‘.ylinder is spread by diffraction (Lind, 1966, 


Greenberg, 1974). In this case eq. 3.2 was used to compute m . The shape of 
the angular distribution is quite similar. When the scattering tensor entering 


in eq. 3.9 was replaced by the exact solution for infinite cylinders normalized 


to "the Rayleigh approximation for spheroids (eq. 4.6), the solution using the 
Rayleigh approximation alone was closely reproduced but the experimental values 
were further approached. This was expected as a comparison between the shape 
of the scattering from the cylinder to the Rayleigh approximation (Fig. 4.7) 
illustrates the well-known fact that back-scattering is overestimated in the 
Rayleigh approximation. 
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: at various scattering augles 8. 


I (linear scale, arbitrary units) 






Fig. 4.2b and c 


Scattering intensity versus azimuthal angle for an ensemble 
of randomly aligned x « 0.47 acrylic cylinders (m “ 1.61) 
with elongation 2:1. The cylinders are distributed vithin a 
"sphere" of radius 5.0 cm which is rotated about different 


axes, 
b) J 

l 


intensity for incident and received radiation polarized 
perpendicular to the scattering plane, 
intensity for incident and received radiation polarized 
parallel to the scattering plane. 






4.1.2 Th^secontJarj^gcnueH^ 

The secondary-scattering approximation derived in see. 3.7 assumes chat the 
near field around a scattering cylinder can he approximated by the far-field 
solution. It is seen from Fig. 4.S that the secondary scattered intensity is 
overestimated by this approwh. Thin ucnclusion is further supported by the 
fact chat wo boo the deep dip in f around 90° that would tend to be filled 
in by multiple scattering. Note that the laboratory has a finite nngulnr 
resolution. 


4,1.3 Remarks and conclusions for T41A. 


The individual cylinders have mirror symmetry In the scattering plane and 
thus do no** by them self give rise to eross-polnrization. The "bird's-nest" 
does not possess such symmetry, but the cross-polarization components were 
too small to be detected in the laboratory (<!). In analogy to extinction 
by a very large object, such as an apple in visual light, the extinction 
cross section in the far field (twice the geometrical cross section for the 
apple) seems to be larger than in the near field (once the geometrical cross 
section). We thus tend to over-estimate the importance of multiple scattering. 


The 


o ! nal n-cna b fr /sv* trier 
“■ * * *© * " Cl 


nnn rav imip (mt tmidnuor 

jr jr **'^*r-.* *'’■*»***■*' fr #«r 




to satisfactory reproduce 


the experiment especially in tin b..(k~scnl taring hemisphere. At small 
scattering angles, Mie-theory (-Fraunhofer diffraction' gradually becomes a 
better approximation. The degree of polarization shown in Fig. 4,8 exhibits 
a number of oscillations at 0 < Sf»° that are attributed to the biref ringent 


character oi the ensemble. 
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0 * 60 ’ 120 ° 180 * 

0 


Fig. 4.3 Theoretical and experimental intensities versus scattering angle 

for the ensemble T4IA of 4:1 elongation aligned cylindrical particles 
with size parameter ?.xa/A ■ 0.47. I is the intensity for radiation 
polarized perpendicular to the scattering plane. Dotted curve is 
the arithmetic sum of individual particle scattering intensities; 
solid curve is the Mie theory result for the ensemble treated as a 
homogeneous sphere with an equivalent index of refraction (see 
Equations 4.7 - 4.9); dashed curve is for the ensemble scattering 
treated coherently and averaged statistically using the single scattering, 
"interference" approximation (see text for explanation) the scattering 
by individual cylinders has been approximated by that of equal volume 
spheriods using the Rayleigh approximation; circles are the experimental 
results. 
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Fig. 4.5 Same as Fig. 4. 3 except in the theoretical computations the 

finite cylinders are approximated by equal length portions of 
infinite cylinders. The dots includes the secondary scattering 
approximation, triangles are the single scattering approximation 
and the solid curve the mie-computations using an index refraction 
given by eq. 3.2. Squares are the experimental results. 
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Fig. 4.6 Same as Fig. 4.5 except I, 2 is the intensity for radiation 
polarized in the scattering plane. 
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Fig. 4.7 
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I O * O 1 < 

120 150 180 


Theoretical intensities versus scattering angle used to approximate 
the scattering by the finite cylinders at perpendicular incidence 
making up the ensemble T41A. For comparison the exact solution for 
the spheriod »7as computed using a computer code provided by R. 
Schaefer for details we refer to R. Schaefer, i960. 
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Fig. 4.8 Theoretical and experimental degree of polarization versus scattering' 
angle for the ensemble T41A composed of 125, 4:1 elongation aligned 
cylinders. Solid curve i6 the arithmetic sum of scattering by 
individual particles; dashed curve is the Mie theory remit for the 
, ensemble treated as a sphere with a corresponding isotropie index 
of refraction obtained through eq. 3.12 dots are for the ensemble 
scattering approximated by the single scattering "interference" theory 
(see text for explanation) the cylinders have been approximated by equal 
volume Rayleigh spheroids; crosses are the experimental results. 
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' I 

Fig. 4.9 Theoretical (solid line) and experimental (dots) normalized standard 
deviation of scattered intensity versus scattering angle from the 
ensemble T41A of 125, x c 0.47 aligned 4:1 elongation cylinders as 
the ensemble is rotated" in azimuth. is the intensity for radiation 
polarized perpendicular to the scattering plane. 



0 30 60 90 120 150 180 

© 


Fig. 4.10 


Same as Fig. 4.9 except I is the intensity for radiation polarized 
in the scattering plane. 
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4.1.4 Scattgring from Pj^l j,_a _two_dimgnslonal ^''bitd’o^ ne st** 

PI Is 156, x ■ 0.47, 1:1 -elongation homogeneous lucite cylinders randomly 
aligned within an X * 19.76 circular disc. 

i 

The scattering o£ each cylinder can be closely approximated by that of a single 
dipole. The simple geometry of the ensemble yield well separated predictions 
by single, multiple and by dependent scattering approximations, i 

Target P1 1 was first orientated such that the plane containing the 
scatterers was perpendicular to the scattering plane. The azimuth angle x refers 
to the angle between the plane of the targets and the incident $ vector (see 
Fig. 4.11). 



ks 

Fig. 4.11 Orientation of target P 1 1 during measurements presented in Figs. 

4.12 and 4.13. The plane containing the scatterers is perpendicular 
to the scattering plane t 0 t s (YZ) , x is the angle in the scattering 
plane between the incident radiation H Q and the plane containing 
the scatterers. 


In Fig. 4.12, the averages I (x) and I (x) over are represented by 

11 22 * 

the solid and dashed curves, respectively, at 9 ® 45° and the same quantities 
at 0 a 90 are shown in Fig* 4.13.^We se^ that in all cases a maximum occures 
at x a At this orientation, k n « k n where u is the normal to the plane 

containing the scatterers. This is the well known condition for reflection 
from an infinite plane. As the plane is finite in our case, we expect inter- 
ference patterns which simply correspond to the diffraction pattern that would 
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be obiBexved if the plane was homogeneous. The scattering pattern is therefore 

symmetrical around the maximum. At 0 * 90°, the individual cylinders have a 

very sharp dip in I . Thus, obs rvations of I (x) at 0 * 90° yield information 

2 2 2 2 

about the importance of multiple and dependent scattering. The x * 0/2 

orientation corresponds to the Brewster angle of an infinite plane. Thus, 

rather than the maximum at this orientation being seen for J , dependent 

scattering predicts a minimum in J . Dependent and multiple scattering are 

also expected to give rise to cross-polarization when the agglomerate does not 

possess mirror symmetry in the scattering plane. This occurs when the plane of 

the scatterers deviates from the scattering plane and does not contain the 

perpendicular to this plane. Significant cross-polarization could not be 

detected at any orientation, and it is seen in Fig. 4,13 that there is no dip 

in the J component at x " 45° and 6 “ 90°. We conclude that multiple and 
22 ... 
dependent scattering are much weaker than single scattering for this target. 





LOGjq (I ) ( I is in arbitrary units } 
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Fig. 4.12 Scattering averaged over ¥ at 0 = 43° by the two dimensional 
ensemble of scatterers P11 shown in Fig. 4.11 as the ensemble 
is rotated in azimuth around the X-axis, solid curve is the 
component of scattered intensity; dashed curve is the I ( 
component, 
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Fig. 4.13 Same as Fig. 4.12 except 8 « 90°. 



4.1.5 Sca ttering, from U41A and C41 

U41A; 43, x • 1.88, 4 ; 1 -elongation homogeneous cylinders aligned parallel 

to each other within an X ■ 25,09 sphere. 

<341; 1, x - 1.88, 4:1 -elongation homogeneous cylinder. 

The differences between T41A and U41A are: 

I) The cylinder size parameter x has been increased a factor 4, such that scatter- 
ing firotn the individual cylinder C41 could be detected in the laboratory. 

II) The "hird's-nest" size parameter X has been increased a factor 2,5. 

III) The average polarizability per unit volume as determined by eq. 3.1, is 
increased by fitting 43 cylinders into the ’’bird's-nest" rather than 30.5 
in the same volume. 

It is evident from Fig. 4.14 that the ensemble does not act like a sum of 

the individual cylinders. The scattered intensity per cylinder is much lower 

for fciie agglomerated cylinders than for independent identical cylinders C41A 

in the range of 0 covered by the experiment. A sharp rise in both J and J ^ 

is expected both from constructive interference and from Fraunhofer diffraction 

at smaller B’b, This target seems to be well suited for investigation of 

multiple and/or dependent scattering effects that may have flattened out the 

angular distribution. As a result, the run of the degree of polarization with 

0 (Fig. 4.15) is flatter and generally lower than for individual cylinders. 

In an attempt to separate an eventual enchancement in J independent of the 

cylinder orientations, U41A was aligned such that the cylinder axes were 

•parallel to the scattering plane (iji » 90°) and x “ 60°, The resulting angular 

distribution of the scattered light (Fig, 4.1 6b) has a double peak in I ^ around 
o 

0 “ 100 whereas a single peak can be seen in I . However, the observed peaks 

22 

are much too weak to fill in the dip in I at perpendicular incidence and 
o 

6 Rt 105 , For comparison the scattering by the similarly orientated single 
cylinder C41 identical to the 43 cylinders in U41A is shown in Fig. 4.16a. We 
note that the peak around 0 « 2x = 120° characteristic of elongated targets, 
has almost vanished for the spherical ensemble. 

In conclusion, it seems probable that the scattering by the ensemble is 
smoothed out by multiple scattering with respect to that by individual cylinders. 
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Fig. 4.14 Intensities per cylinder versus scattering angle for x * 1.88 cylinders. 

Dashed and solid curves are respectively computed J, and components 
for a 4:1 portion of an infinite cylinder at perpendicular incidence. 
For a‘ similarly orientated single 4:1 -elongation cylinder crosses and 
dots are the experimental values for J and I 22 respectively. The 
triangles are the experimental I component for the ensemble U41A, 
circles the I 22 component. 
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Fig. 4.15 Degree of polarization versus scattering angle to..* x “ *.88 cylinders. 
Solid curve is computed ^or an infinite cylinder at perpendicular 
incidence. Experimental results are: pluses » single 4:1 cylinder at 
normal incidence (target C4‘1 A) , circles are for the ensemble U41A of 43 
aligned 4:1 cylinders. 
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0 

Fig* 4.16a Scattered intensities in units of the potential created over 
the detector versus scattering angle, from the 4:1 -elongution 
x - 1,88 cylinder C4 1 The axis of symmetry is contained in the 
scattering plane and makes an angle <p •* 60° with the direction 
of the incident radiation k 0 , During this experiment the back- 
ground was not compensated for. The uncompensated background 
radiation is displayed in the same units in append! * . 



0 

Fig. 4.16b Scattered intensities in units of the potential created over 
the detector versus scattering angle, frem the 43 aligned x ** 
1.88, 4:1 -elongation cylinders confined into the ensemble U41A. 
ip m 60° for each individual cylinder. During the experiment the 
background was not compensated for. The uncompensated background 
radiation is displayed {.n the same units in appendix E. 
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4.1.6 U^R_and > _CAjtR 


U41R: 43, x ■ 1.80, 4.* 1 -elongation homogeneous cylinders randomly aligned 

within an X * 20.08 sphere, 

C41R: 1, x ■ 1.88, 4:1 -elongation homogeneous cylinder, the averages are 

taken over random orientations, 

The differences between U41R and U41A are: 


X) The 43 cylinders are not only randomly located as in U41A but also randomly 
orientated rather than aligned ’perpendicular to the scattering plane, 

II) The size of the spherical ensemble is decreased from X « 25.09 to 20,08, 

The difference between C41R and C41A is: 

I) The measurements C41R are performed over uniformly distributed orientations 
to simulated randomness. The actual target C41 is identical in the two 
experiments, 


The angular dependence of the two polarization components I and J in 

XX 2 2 q 

experiment C41R are quite similar and approximately straight lines up to 0 * 90 

on the logarithmic plot in Fig, 4,17. Althr jgh very shallow, the dip at 0 w 105° 

in I can still be seen. Otherwise the scattering curve is relatively flat with 
li 0 

very little backscattering up to the last measurement ?> 8 ■ 165 , This 
scattering curve is discussed in more detail by Schuerman at al, (1980) and compared 
to those for a selection of 28 targets of various size and shape. As the cylinders 
in target U41R are randomly aligned, the resulting averaged index of refraction 
over a portion of the agglomerate is independent of the polarization of the 
incident and scattered light. Thus, the single-scattering approximation predict,'? 
the Bame changes in all polarization components (shown in Fig. 4.18), Higher- 
-order scattering and/or dependent scattering would add a polarization-dependent 
phase shift with respect to the single scattering. The difference in the degree 
of .polarization of light scattered from the single cylinder and from the agglomerate 
respectively Fig. 4.19 thus illustrates the importance of these effects. 
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Fig. 4.18 Measured scattered intensities per cylinder versus scattering angle 
for 43, x = 1.88, 4:1 -elongation cylinders randomly located and 
orientated in the X **■ 20.08 spherical ensemble U41R averaged over 
uniform orientations, of the agglomerate. 
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Fig. 4.19 Measured degree of polarization averaged over uniform orientations 
versus scattering angle for the x = 1.88, 4:1 -elongation single 
cylinder C41R and for the ’'spherical" X = 20.08 ensemble U41R of 
43 such cylinders. 
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4.1.7 Dep endence on cylinder elongation 


T41R: 125, x ■ 0.47, 4:1 -elongation homogeneous cylinders randomly aligned j 

within an X ■ 10.04 sphere, 1 

•B41R: 125, x * 1.74, 4:1 -elongation homogeneous cylinders randomly aligned 

within an X ■ 37.25 sphere, 

T21R: 240, x ■ 0.47, 2:1 -elongation homogeneous cylinders randomly aligned 

within an X * 10.04 sphere. 

B21R: 240, x <* 1.74, 2:1 -elongation homogeneous cylinders randomly aligned 

within an X » 37.25 sphere. 

T11R: 539, x » 0.47, 1:1 -elongation homogeneous cylinders randomly aligned 

within an X ■ 10.04 sphere. 

B11R: 539, x ■ 1.74, 1:1 -elongation homogeneous cylinders randomly aligned 

within an X * 37.25 sphere. 

T41R, T21R and T11R represent three agglomerates differing in the length- 
-to-diameter ratio of the individual cylinders, the average polarizability is 
kept constant by varying the number of cylinders. B41R, B21R and B 1 1 R are the 
same targets, respectively, but the measurements were performed at the 
millimeter-wave laboratory at Bochum. T41R and B41R are compared in the next 
section (4.1.8) for colour effects. It seems from Figures 4.20, 4.21, 4.22 and 

4.23 showing the scattering functions as a function of 0 for T41R, T21R, T11R 
and B41R, B21R and B11R, respectively, that the cylinder elongation has a 
negligible effect on the scattered intensities averaged over all target orienta- 
tions. However, the dependence on orientation described by the standard deviation 
o of the 272 measurements should be a measure of the target isotropy. In Fig. 

4.24 the normalized standard deviation a/I is plotted as a function of 9. It is 
seen that a/I generally decreases with the length of the silicate rods, which 
are thus more evenly distributed. Scattering in the forward direction is 
dominated by the Fraunhofer diffraction, which is almost independent of the 
orientation of the spherical agglomerate. We note also that the phase of the 
scattered radiation from individual cylinders and therefore the resulting inter- 
ference pattern depends most heavily on m and the size of the agglomerate at 
sufficiently small scattering angles. In the case discussed here. X(/n - 1) 

dC X, and therefore the interference is generally constructive; i.e., the 
dependence on target orientation is small. At larger scattering angles, on the 
othpr hand, the phase depends more heavily on locations of the particles within 
the ensemble. It is therefore widely differing from particle to particle. The 
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interference pattern as the target T21P. is successively rotated around a few 
axes perpendicular to the scattering plane is shown at 135° scattering angle 
in Figures 4* 2b and 4,2c* on p. 42. 
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Fig. 4.20 Measured scattered intensities averaged over uniform orientations 
versus scattering angle for three X = 10.04 "spherical" ensembles 
of randomly orientated x = 0.47 cylinders. There are 125, 4:1 
-elongation cylinders in target T41R, 240, 2:1 -elongation cylinders 
in T21R and 539> 1:1 -elongation cylinders in T11R. I were too 
small OS 0.1) to give reliable data. 12 







Pip 4.21a Measured scattered intensities averaged over rotation in azimuth 
about five different axes of the target in experiment BUR. The X 
- 37.25 ensemble of 539 , x - 1 . 74, 1*1 -elongation randomly located 
and orientated cylinders is the same target as T11R but irradiated 
by the X * 8.6 mm (35 GHz) waves of the p-wave laboratory at Bochum 
rather than the 3.18 cm (9.4 GHz) at Albany. When dots and crosses 
coincide only crosses are shown, similarly only squares are shown 
when they coincide with circles. 
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Fig. 4.21b Measured run of the degree of polarization in experiment B11R. 
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pig, 4.22a Measured scattered intensities averaged over rotation in azimuth 
about five different axes of the target in experiment B21R. The 
X - 37.25 ensemble of 240, x « 1.74, 2:1 -elongation randomly 
located and orientated cylinders is the same as target T21R but 
irradiated by the X * 8.6 mm (35 GHz) waves of the p-wave laboratory 
at Bochum rather than the 3.18 cm (9.4 GHz) waves at Albany, when 
dots and crosses coincide only crosses are shown, similarly only 
squares are shown when they coincide with circles. 
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Fig. 4.22b Measured run of the degree of polarization in experiment B21R. 
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Fig, 4,23a Measured scattered intensities averaged over rotation in azimuth 
about five different axes of the target in experiment B4‘1R. The 
X » 37.25 ensemble of 125, x * 1.74, 4:1 -elongation randomly 
located and orientated cylinders is the same target as T41R but 
irradiated by the X « 8.6 mm (35 GHz) waves of the u-wave laboratory 
\ at Bochum rather than the 3.18 cm (9.4 GHz) waves at Albany. When 

3 dots and crosses coincide only crosses are shown, similarly only 

squares are shown when they coincide with circles, 
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Fig. 4.23b Measured run of the degree of polarization in experiment B41R. 
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4.1.8 Cg_^ |jg__de pa nd onc6 i f or, tlic agglqmQcabes of Bnicntea 


BAIR: 125, x - 1.74, 4:1 -elongation homogeneous cylinders randomly aligned 

within an X » 37,25 sphere. 

U41R; 43, x * 1,88, 4:1 -elongation homogesieous cylinders randomly aligned 

within an X * 20.08 sphere, 

T41R: 125, x ** 0.47, 4:1 -elongation homogeneous cylinders randomly aligned 

within on X » 10.04 sphere, 


Experiments T41R and B41R were performed using the same target but at 

different wavelengths. Thus, the differences arc, save for a differing averaging 

technique, truly colour effects. The measurement B41R was performed at the 

laboratory at Bochum operating at 35 GHz (8.57 mm), If we let T41U represent 

a particle in yellow light, B41R would represent the same particle near the 

wavelength of 'Lyman a. The dip in I at the longer wavelength of T41R in Fig. 

? 2 

4.20 is only reduced with respect to chat of single cylinders by multiple or 
dependent scattering. At the shorter wavelength of B41R, however, the single 
cylinders exhibit only a very shallow dip which vanishes in the scattering of 
the "bird's-nest" Fig. 4.23, The phase function for x ■ 1.74 and x * 1*88 
cylinders are not expected to differ very much. The total brightness for 
infinite and randomly orientated cylinders is show* in Fig. 4.25 and for the 
agglomerates in Fig, 4.26. Thus a comparison between the B41R and U41R measure- 
ments Figures 4.23 and 4.18 merely reflect the difference in agglomerate size 
and a small change in the average index of refraction as 43 rather -than 30.5 
cylinders were fitted in the X * 20.08 hollow eccofoam sphere to make up 

"bird's-nest" U41R, There are few differences in the runs of I (0) and I (0) 

11 22 

as given by the two sets of measurements, but the relative importance of cross- 
-polarization 'in B41R is surprisingly low. No measurements for the x « 1.74 
single cylinder was performed to provide an adequate comparison, note also that 
the cross-polarized components were not rigorously calibrated as these 
components vanish for the spheres used as standard (see also section 2.5.6). 


We conclude that in this case the size of the agglomerate docs not seem 
to have a decisive effect on the shape of the scattering function, and that the 
magnitude is approximately proportional to the number of scattering particles. 
Whereas the colour dependence of the individual scatterers is still quite 
important for agglomerates of this size the effect is expected to be reduced by 
dependent scattering which increases in importance as the size of the agglomerate 
grows . 


arbitrary units) 
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Fig. 4,26 Measured angular distribution of brightness per cylinder for the 

ensembles T41R consisting of 125, 4:1 -elongation x “ 0.4? cylinders 
in an X ■ 10.04 "sphere", 04 1R; 43 4:1 -elongation x » 1.88 cylinders 
in an X “ 20,08 "sphere" and B41R consisting of 125 4:1 -elongation 
x => 1,74 cylinders in an X * 37.25 "sphere" respectively. 
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4. 1 .9 (Depandence on particle number density, 

T21R: 240, x « 0.47, 2:1 -elongation homogeneous cylinders randomly aligned 

within an X * 10,04 sphere. 

1 TC21R: 204, x » 0,47, 2:1 -elongation homogeneous cylinders randomly aligned 

within an X «■ 15,81 sphere. 

V21R: 44, x ■ 0,78, 2:1 -elongation homogeneous cylinders randomly aligned 

within an X ■ 15,81 Bphcre. 

The Fraunhofer diffraction concentrates towards the forward curection as 
the size of the ensemble Increases, It is thus expected that the side scattering 
from TC21R is more closely represented by the pure "interference" approximation 
than the more' closely*' packed T21R, This approximation also predicts a concentration 
towards small 0'a for the expanded ensemble (see section 3,6,3). The dependence 
on orientation x» <1*, (a, 6) and therefore a/I is expected to be smaller for TC21R 
at GSO but to iwarease more rapidly with 0 (fig. 4.28) whereas the resulting 
scattered intensity (fig, 4,27) decreases more rapidly. This can bo seen from eq. 

3.21 in the limit m^-M Target V21R occupies the same volume as TC21R, but the 
scattering material has been concentrated in larger cylinders. The side scattering 
pattern, (fig, 4.27), mostly reflects the difference in scattering functions for 
the individual cylinders, This can be seen by comparing the total brightness 
measured for TC21R and V21R (fig, 5.5) to the computed brightness for the respective 
cylinders (portions of infinite cylinders) independently scattering the incident 
light (see fig. 4.25), The "fluctuations" in the phase of the scattered signal from 
each cylinder as a function of orientations is expected to vary more heavily for 
V21R than for the more "homogeneous" TG21R. That effect is seen in Fig. 4.28 as 
a flatter 0 dependence for V2IR than for TC21R, 












0 (degrees) 



Pig. 4.28 Measured and normalized standard deviation of scattered intensities 
versus scattering angle for targets; T21R consisting of 240 x » 0.47 
2:1 elongation cylinders, TC21R:204 x = 0.47 2:1, V21R: 44 x = 0.78 
2:1 elongation cylinders in X ® ?'.'.04, 15.81 and 15.81 "spherical" 
ensembles respectively. 
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4.2 Silicate Core - "Ice" Mantle Cylinders 

The scattering by the ensemble IC3b, is compared to that by lC4b and XC4c 
for colour .effects in sections 5.2.2 and 5,2.3, The ensemble ICIa is also 
discussed in the same sections, The run of the polarized intensity 
components are shown separately in 'figures 4.37 and 4.38 for IC3b and 
ICIa respectively. 

4.2.1 PgEgndence_on_the , .size of ...the.eggloinQratfe 

XC4d: 80, x - 2.45 (a » ,07), 1.96 (cr r . « .04) to 1 mean elongation eccofoam 

1.7 mantles and x t » 0.78 ~ 2 to 1 elongation luclte core cylinders 

0 

randomly aligned within an X « 19.3 sphere. 

lC4c: 44, x « 2.45 (o m .06), 1,97 (o * .04) to 1 mean elongation eccofoam 

w m x c 

1.7 mantles and x c « 0.78 2 to 1 elongation lucite core cylinders 

randomly aligned within an X ® 15.8 sphere. 

lC4b: 23, x^ « 2.45 (o x - .06), 1,96 (o >* .05) to 1 mean elongation eccofoam 

1.7 mantles and x fi » 0.78 ~ 2 to 1 elongation lucite core cylinders 
randomly aligned within an X 05 12.6 sphere. 

12Bs 1 , x^ » 2.44,1.954 to 1 elongation eccofoam 1.7 mantle and a x c ■= 0.78 ~ 2 
elongation lucite core cylinders. The averages are taken over random 
orientations. 

The run of the scattered polarized components X . X and I with scattering 

X 1, 2 P X 2 

angle was measured for the single core-mantle cylinder "12" representing a 
single interstellar grain. The scattering averaged over "uniformly" distributed 
orientations (as described in section 2.7.1) to simulate randomly orientated 
cylinders is shown in Figure 4.29. The single scattering approximation described 
in chapter 3 is applied to the equal density agglomerates of 23, 44 and 80 grains 
respectively, all of which are reproductions of "12", The average scattering 
efficiencies over random orientations were substituted for that of each of the 
randomly orientated cylinders. The detected signal from "12" at some combinations 
of polarization, scattering angle and orientation were close to the noise level 
( « lO” 2 ) although version C of the laboratory apparatus was used. The thus 
introduced errors propagate to the single scattering approximation (see 
Chapter 3) . It is possible that the oscilliations seen in the computed run 
of I around 110° scattering angle are caused by this effect. The dots in 
Figures 4.30 to 4.35 are the results of the single scattering computations 
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as well as the open triangles, circles and squares in Figure 4.36 the dashed 
curves corresponding to the dots are added for illustration purposes only, as 
are the curves in Figure 4.36, It is seen in the Figures 4.30 to 4.35 that the 
scattering at small angles approaches the diffraction pattern discussed in section 
3.3. The diffraction pattern is approximated by that for a homogeneous sphere 
with "equivalent" index of refraction as given by Equation 3,1 after averaging 
the dipole moments per unit volume (P) . Towards higher scattering angles the 

f 

single scattering approximation is approached. The cross-polarization components 

J * I over uniform distribution of orientations are compared in Figure 4.36, 

12 21 

the diffraction does not polarize, It is therefore expected that the experimental 
results follow the single scattering approximation at all scattering angles. 

Note that the measured scattering by the single cylinder "12", on which the 
"interference" theory relies, might be subject to large experimental errors. 

It was concluded in section 4 . 1 . 1 that the scattering by the 125 aligned 
cylinders in target T41A were closely reproduced by the single scattering 
approximation at 6 > 90° and approach the diffraction pattern at lower 0's. 

We may draw the same general conclusion for targets IC4b, IC4e and IC4d 
except that the magnitude of the I component is lower for the larger 
agglomerates than predicted by the theory. This leads to a lower degree 
polarization for agglomerates than for individual cylinders and moves the 
maximum degree of polarisation do lower scattering angles for the bigger 
agglomerates . 


LOG.n (I) 


Measured intensities scattered off a single silicate-core- 
"ice ,, -mant.;le cylinder of 2:1 elongation modelling an inter- 
stellar dust grain. The intensities are averaged over 
"uniformly"' distributed orientations {see section 2.7.1 for 
explanations). The connecting curves are added for illustration 
purposes only. 
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Fig. 4.30 Theoretical and experimental intensities versus scattering angle for 
ensemble IC4b consisting of 23 randomly located and orientated 
core-mantle cylinders with x = 0.78, x = 2.45 and close to 2:1 
elongation within an X = 12.6 c "sphere". is the. intensity for 
radiation polarized perpendicular to the scattering plane. Solid . 
curve is the Mie theory result for the ensemble treated as a homo- 
geneous sphere with an equivalent index of refraction (see equation 
3.1); dashed curve is for the ensemble scattering treated coherently 
and averaged statistically using the single scattering "interference" 
theory (see text for explanation) the scattering by individual cylinders 
was determined experimentally and is shown separatly in Fig. 4.29; 
crosses are the experimental results for the ensemble. 
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Fig . 4.31 Same as Fig. 4.30 except: I is the intensity for radiation polarized 
in the scattering plane. 
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Theoretical and experimental intensities versus scattering angle for 
ensemble IC4c consisting of 44 randomly located and orientated core- 
-mantle cylinders with x «* 0.78, x *» 2.45 and close to 2:1 elongation 
within an X'= 15.81 "sphere". I is the intensity for radiation 
polarized perpendicular to the scattering plane. Solid curve is the 
Mie theory result for the ensemble treated as a 'homogeneous sphere 
with an equivalent index of refraction (see equation 3.1); dashed 
curve is for the ensemble scattering treated coherently and averaged 
statistically using the single scattering "interference" theory (see 
text for explanation) the scattering by individual cylinders was 
determined experimentally and is shown separately in Fig. 4.29; crosses 
are the experimental results for the ensemble. 
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Fig. 4.1.4 Theoretical and experimental intensities versus scattering angle for 
the ensemble XC4d consisting of 80 randomly located and orientated 
core-mantle cylinders with x„ = 0.78, x « 2.45 and close to 2*1 
elongation within an X «* 19.29 "sphere" 1 '. 1 1 is the intensity >df 
radiation polarized perpendicular to the scattering plane. Solid curve 
is the Mie theory result for the ensemble treated as a homogeneous 
sphere with an equivalent index of refraction (see equation 3.1); 
dashed curve is for the ensemble scattering treated coherently and 
averaged statistically using the single scattering "interference" 
theory (see text for explanation) the scattering by individual 
cylinders was determined experimentally and is shown separately in 
Fig. 4.29; crosses are the experimental results for the ensemble. 



0 ° 30 ° 60 ° 90 ° 120 ° 150 ° 180 ° 

0 

Fig. 4.35 Same a3 Fig. 4.34 except is the intensity for radiation polarized 
in the scattering plane. 
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Fig. 4.36 Theoretical and experimental results for the cross-polarized intensities 

versus scattering angle. The open triangles are for the ensemble scattering 
treated coherently and averaged statistically using the single scattering 
"interference" theory (see text for explanation). 



ig. 4.38 Measured polarized intensity components ve.sus scattering angle for 
target ICta consisting (i 243 x = 0.47, 2:1 elongation lucite 
cylinders and "waste" material from the manufacturing of the eccofoam 
1.7 (m fa 1 .3) mantles, in an X = 10.04 "spherical" ensemble. 
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Chapter 5 


COMPARISON WITH ASRTRONOMICAL OBSERVATIONS 


5.1 Introduction 

Observations of comets as well as of the zodiacal light yield information 
about the sura of scattered intensities off the contributing scatterers, integrated 
along the line of sight. The geometry of the system in which the observations 
are made are related to the laboratory system as shown in fig. 5.1, The translation 
of observable quantities between the two systems is not readily made and is 
discussed elsewhere (see Dumont 1976 for references), e is the angle between the 
direction to the light source and the line of sight. The plane defined by the 
location of the observer, scatterer and the Sun which is approximately a point 
source, is the scattering plane with respect to which the direction of polarization 
is defined (see section 2.3). 


cnti 

sun 



Fig. 5.1 ' The geometry for observations of scattering by particles 

in the interplanetary medium. 
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The widest coverage in angular (6) distribution is found in zodiacal light 
measurements. Translation of the observations into quantities that can be 
measured in the laboratory involves assumptions regarding the particle number 
density distribution in the scattering plane. It has traditionally and for 
simplicity been assumed that the particle number-density (») is proportional 
to a power of heliocentric distance 7?(n <r J? -V ) and that the light scattering 
functions are independent of location. As pointed out by Schuerman (1980) r this 
last assumption has no theoretical support whatsoever. On the contrary, a 
segregation is expected as a function of heliocentric distance. Bearing in mind 
that these assumptions are not verified by observations from the Pioneer 10 
Spacecraft (Schuerman, 1980, refered to above) we might compare the angular 
dependence of the brightnass (J ) and degree of polarization (P) of the 
zodiacal light per unit volume to the laboratory data for our models of cometary 
debris. The most complete surveys of the zodiacal light along the ecliptic 
were transformed to the laboratory geometry by Dumont (1976). 

5.2 Comparison to the Zodiacal light 

The zodiacal light has been found to have approximately the colour of the Sun 
(for references see Weinberg and Sparrow, 1978). We are therefore first discussing 
the general characteristics of the zodiacal light in sections 5.2.1 and 5.2.2, 
while colour effects are discussed in section 5.2.3. Figures 5.2 and 5.8 are 
reproduced from Dumont (1976) by permission from Springer-Verlag. 

5.2.1 Brightness 

The total brightness (I * + I lz + -P 22 + J 2l ) deduced by Dumont from measure- 

ments at Haleakala (Weinberg, 1964) and Tenerife (Dumont 1965, Dumont and Sanchez 

1975), is shown in Fig, 5.2 where v ** 1.? was used. I. ^ reaches a local maximum 

tot 

at fl « 33° and* becomes roughly isotropic from 50° to 130° after whi«h it rises in 
backscattering. The local maximum at 33° scattering angle is shown in greater 
detail by Dumont (1976 fig. 5). 

Although appearing in both the Tenerife and the Haleakala data, and making 
a smooth connection to the rocket data by Leinert et al. (1974) that probes the 
local minimum around 6 * 20°, the very existence of the peak is dubious. This 
is both because the accuracy of the inversion method and the reliability of the 
observations decrease with scattering angle. This region is labeled poorly 
reliable by Dumont. It is, however, noteworthy that a similar peak is seen in 
the laboratory data for target T41A. It is also predicted by the "interference" 
theory for most "lird’s-nest" particles, but the angular resolution of the 
laboratory data i. too low to more than hint at its existence (see figures 5.3 t.o 
5.6). 
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(A second peak around 6»55° in the data for TftlA is also seen in Dumont's 
Fig. 3, 1975.) Dumont (1976) writes; 

"The variations of a from 8 « 50° to 0 * 130° are probably within 
a factor of 2, and this conflicts with many theoretical scattering 
functions which show variations by a factor of 5 to 10 in the quoted 
range of scattering angle." 

a is the total brightness (J in this text) . Most of the laboratory data were 
taken with 15° intervals, so that J vas measured at 6 = 45°, 60° and 120°, 135° 
but not at 50° and 130°. It is seen in table 2 that scattering by "bird' s-nests" 
is generally closer to isotropy in this range than scattering by individual cylinders 
and that the agglomerates of cylinders with eroded mantles have a flatter run 

of X ' . with 9 in this region than the agglomerates of "ice” coated cylinders. 

tot o o 

Dumont estimates the variation in the whole interval of 0 s 20 tu 180 to be 

within a factor of 5, and writes; 

"Such a flatness can agree with few of the scattering functions repor- 
ted by Wickramasinghe (1973) for Mie particles, or by Giese (1970, 1971), 

Giese et. al. (1974) for elaborate a.ixtures of homogeneous or mantle-core 
particles: most of these curves exhibit a much stronger ft_*ward scattering. 

Some agreement may perhaps be sought with absorbing particles, for which 
the ratio a la. is generally of the order of 10 in the 20°< 0 < 180° 
range, but can be as low as 2 or 3. On the other hand, the curves of figs 
4 and 5 seem very hard to reconcile with dielectric particles, for which 
the same ratio is between 20 and 200 or more." 

Fig. ft in Dumont's article is here Fig. 5.2, the 10° to 60° range is shown in more 
detail in Dumont's Fig. 5. Although none of the scattering functions for the 
"bird 1 s-nests" actually investigated in the laboratory show quite such an isotropy 
in this range of scattering angle, we note that some agglomerates that are dielectric 
as are all "bird -nests" reported on in this study, show J /I . <5. 

Agglomerates of more realistic sizes for typical zodiacal-light particles are 
currently being studied, using the "interference" theory. 

In the backs c&ittering region 160°< 0 <; 180° there is a rise in J (Fig. 5.2). 

Dumont estimates the ratio I^ rt ^(180 °)/ 1 . to 1.7 for the Tenerife data 

tot tot, mm 

and to t.9 for the Haleakala data. A rise with a corresponding slope is evident 
in several of the scattering functions for "bird' s-nests", where 1 < ■^ tot (‘l65°) / 
u tot miir^ ^ (see Fxgs. 5.3 to 5.6) . 

In agreement with Dumont's Fig. 6 in the same reference, the component 
( X in Dumont's Fig.) is flatter than the I 2 component (//) for all "bird 1 s-nests" 
in this study. 


type of cylinders! 

Target Id. § single agglomerate homogeneous core-mantle rang* in 8 T / J ^ . 

(decrees) tot ’ max 

less than 
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Fig. 5.2 Reproduced by permission given by Springer-Verlag, from Dumont 
(1976, Fig. 4, p. 92). The phase function for interplanetary 
dust deduced from the Tenerife (Dumont 1965, Dumont and Sanchez 
1975) and Haleakala (Weinberg, __1 964) photometric surveys. The 
number density n given by n<*R v where R is the heliocentric 
distance and v =» 1.2 was assumed. 
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Pig. 5.3 Measured total brightness per cylinder versus scattering angle 
fdi the single core'-mantle grain M 12 H and the agglomerates IC4b, 
IC4c and IC4d of similar grains. For target parameters see the 
fold out page. The connecting curves are added for illustration 
purposes only. 



e 

Fig. 5.4 Measured total brightness versus scattering angle for 
agglomerates TC21R and T21R differing in the silicate 
grain number-density and for the agglomerate ICIa 
containing a mixture of silicate cylinders and H ice”- 
type material. See the fold out page for target para- 
meters. 
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Fig. 5,5 Measured total brightness versus scattering angle for 
the equal volume "bird* s-nests" TC21H; triangles, V21R; 
squares and IC4c; open boxes. See the fold out page for 
target parameters. 
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Fig. 5.6 Measured total brightness versus scattering angle for 
the agglomerates IC4b, IC3b and IC4c consisting of 
core-mantle cylinders. For target parameters see the 
fold out page. 
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Fig. 5.7 Measured total brightness versus scattering angle for 

the ensembles T41R, T21R and T11R consisting of x-0.47 
silicate type material cylinders, differing only in 
elongation. For target parameters see the flod out page 
The solid line is the theoretical computation of the 
scattered brightness from an equal mass portion of an 
x = 0.47, m = 1.618 - 0.005i, infinite cylinder averaged 
over ’’random" orientations. 
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5,2.2 Degree of Polarization 

The similarly deduced degree of polarization (P(0)) is less sensitive 

i to the value assigned to v than the photometric curves (Dumont and Sanchez , 

1975, fig. 5). In fig. 5.8 P(8)is reproduced from Dumont (1976, fig. 7), 

where the same value forv (1.2) is used as for inversion of the intensities. 

P is generally positive at 8 > 60°, the apparent large negative degree of 

polarization from the Haleakala data in the poorly reliable left part of 

the figure is not believed to be genuine by Dumont. From Dumont's fig. 6 

it appears that both curves join in a neutral point at flw 35°. The maximum 

degree of polarization, P roax & 0.4 occurs close to 0 «* 90° in the Haleakala 

inversions and closer to 80° for the Tenerife data, with P max w 0.35. The 

values of these curves at 6 « 90° and 0 = e are of special importance 

m 

because they are independent of. v (Dumont 1972, 1973). is the elongation 
of the maximum degree of polarization, 55° - 75°. 

It can be seen in fig. 5,9 to 5.13 that the general shape of these curves 

are reproduced by the "bird’s-nests". P is generally positive, the magnitude 

is close to that measured for the agglomerates IC4b , IC4c and IC4d (fig. 5.9) 

of "ice"~coated silicate, for which the two points P(e OT ) and P(90°) are in 

excellent aggreement with the observations. The maximum however occurs at 

higher 0’s but is expected to shift towards smaller 0's with the Fraunhofer 

di£ fraction" as the ."bird* s-nests” size increases. P Biax for all the agglomerates 

with eroded mantles is substantially higher than those of fig. 5.8 but the 

maximum occurs closer to P in the inverted observations. 

max 
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Fig. 5.8 Reproduced by permission given by Springer-Verlag , from Dumont 
(1976, Fig. 7, p. 95). The degree of polarization versus 
• scattering angle for interplanetary dust deduced from the Tenerife 

(Dumont 1965, Dumont and Sanchez 1975) and Haleakala (Weinberg, 
1964) photometric surveys. A number density n given by n a R where 
R is the heliocentric distance in A.U. and v ■ 1.2 was assumed. 

The left side of the figure is labelled poorly reliable by' Dumont. 
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Fig. 5.10 Measured degree of polarization versus scattering angle for 
agglomerates TC21R and T21R differing in the silicate grain 
number density and for the agglomerate ICIa containing a mixture 
of silicate cylinders and "ice ,r -type material. See the fold out 
page for target parameters. 
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Fig. 5.12 Measured degree of polarization versus scattering angle for the 

agglomerates ICAc, X,C3b and IC4b, of ''core-mantle" cylinders. For 
target parameters, see the fold out page. 
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Measured degree of polarization versus scattering angle for 
the ensembles T41R, T21R and T11R of x = 0*47 silicate 
cylinders, differing only in elongation, for target para- 
meters, see the fold out page. • ' ' 
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Fig. 5.14 

•Measured degree of polarization versus scattering angle for 
single silicate-core - "ice"-mantle cylindrical model of an 
stellar grain. See the fold out page for target parameters. 

the 

inter- 


5.2,3 Co lour 


If Che diameter 0,1 Mm is assigned Co Che silicate cores, targets lC4b, 

IC4c, IC4d and V21R models scattering of X " 0.4 ym light, IC3b scatter at 
X - 0,5 ym and ICIa, T11R, T21R, T41R and TC21R all scatter 0.7 ym light. 

.Pure colour is represented in fig. 5,6, both IC4c and IC3b correspond to 2 ym 
diameter agglomerates of 44 interstellar grains. The colour is generally blue, 

I fl » 3. If the agglomerates reach sizes that are currently believed 

0 H 0 • 5 

to* be more realistic for typical zodiacal-light particles (10-100 ym), the size 

parameter of the individual particles might be optically more important than the 

size parameter of .he agglomerates* Note that in the limit where the particles 

become homogeneous, scattering by Mie-particles with a given index of refraction 

is almost colourless in this size range. We have therefore also compared the 

scattering from ensembles of equal-size parameter X while the size x of the 

individual grains is varied. J tQt .from IC3b ,and IC4b, both with ensemble size 

parameter X * 12,65 are also compared in fig. 5,6 where it is seen that the 23 

larger cylinders scatter more light than the 44 smaller grains in IC3b, It appears 

therefore that the blue light (X * 0,4 ym) is scattered, more efficiently than 

green (X “ 0.5 ym) light, I /I w 2. The degree of polarization (fig. 5.12) 

0.4 0*5 

is lowest in the blue light (IC4c and IC4b) the position of the maximum 
0 « . w 12u' is shifted to ra 90^ at X ® 0.5 ytn« The agglomerates .TG21R and 

Of 

V21R of bare cylinders both have X * 15.81 in X « 0,7 ym (red) light and X 
- 0.4 ym (blue) light respectively. Similarly to scattering by IC4b compared to 
IC3b,2\ . is higher for V21R than for TC21R, I /J £ 2, suggesting s.hat the 

tOt 0,4 0 # 7 

light scattered off the "ice*' - free agglomerates also is slightly blue (fig. 5,5), 
On the other hand, the run of the degree of polarization looks quite similar for 
both ensembles (colourless) with P mx 0*87 around 0 ■ 90° (fig. 5.11), 


Few differences in, the zodiacal light observed at different visual wave- 
lenghts are larger than variations in the same part of the, spectrum between 

t 

observations by different astronomers. The colour of the zodiacal light 
is .therefore believed to be white or grey. The colours indicated by the 
microwave measurements might be too prominent to compare favourably with the 
observations. As a realistic comparison involves scattering functions over 
a size distribution, more microwave measurements are desirable. It seems, 
however, that the scattering efficiency per cylinder decreases as the 
agglomerates of core-mantle grains grow in size (fig, 5.3). We see in 
figures 4.30 to 4.35 that the T component decreases faster than I ^ resulting 
in a lower degree of polarization, (fig. 5.9), approaching that of the zodiacal 
light. 
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Some astronomers have observed negative polarization at high scattering 
angles (Frey at al,, 1974, Wolstcncroft and Rose, 1967). Weinberg and Mann (1968) 
found that the neutral point near 160° tends to smaller elongations with increasing 
wavelength. Neagative polarization is seen at 0 ** 165° in several of the experimental 
graphs for "bird's-nests”, 


5.3 Comets 

The inversion of the brightness intergral is much simpler for comet data 
than for the zodiacal light because all light observed in a given direction 
(c) comes from acatterers distributed in a narrow interval of the line of 
sight, the scattering angle (9) is essentially constant. To the knowledge 
of the author, only observations of the dust tail of comet Ikeya-Saki, 1965 
VIII, (Weinberg and Beeson, 1976, Matyagin et al., 1968), provide information 
on I and P over an appreciable interval of fl. The angular dependence was 
deduced by scanning the dust tail with a photometer. It follows that in 
addition to the dependence on possible variations in scattering functions 
for particles as a function of location in the tail, affecting P(0), the 
measured brightness J(0) alao depends on the particle number density distribution, 
Weinberg and Beeson conclude from model calculations of the degrees of polarization 
using a Mie-program, that the observations arc consistent with grains in a very 
narrow size-distribution. Although the observed polarization is very different 
from that of light scattered off the single core-mantle cylinder 12 (Fig. 5.14), 
the grain parameters are suggestive of those used for the single core-mantle 
grains in this study. 


Some dust properties can also be inferred from dynamical considerations. 
Moat dust tails can be resolved in streamers or synchrones and in syndynes. 
Streamers coincide with the loci expected for particles ejected simultaneously 
from the nucleus and subjected to a particle-dependent heliocentric repulsive 
force - - — — — m * counteracting the gravitation , (/? is the heliocentric 

Z. * jjCalXIX 

distance ana w g ra £ n is the mass of the grain) . Syndynes correspond to loci for 
continuously ejected particles with % given y, Conclusions on the optical 
properties of the particles can be drown if y - 8 is assumed, where 0 is the 
acceleration resulting from radiation pressure in units of the gravitational 
attraction. 


Following the -model suggested by Whipple (1951) and formulated by Finson 
and Frobstein (1968), Sekanina and Farrell (1980) interpret the structure in 
the dust tu * of comet West 1976 VI as evidence for fragmentation of strongly 
non-spherical particles, The parent bodies resulting from a single burst are 
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lined up along a at raamar, As the particles burst, the fragments line up in etrlAe. 

Sekanina and Farrell (1980) conclude that the repulsive acceleration p was 
slightly lower for the parent bodies than the "average 1 ' on fragments ranging 
from 0,6 to 2,7 times the solar gravitational attraction. With the assumption 
U " 6» this relation, which is difficult to account for with Mie-theory, leads 
Sekanina and Farrell (1980) to suggest that the parent particles are chain-liki 
agglomerates* We note, however, that 8 is generally not acting in the radial » 

direction if the scatterers docs not posses spherical symmetry as for elongated 
structures. This is easily realised by noteing that the scattering pattern is 
asymmetric with respect to the incident ^-vector for particles which are 
asymmetric with respect to this vector. We note also that the formation of 
discrete striae requires a distinct correlation between the fragmentation 
media »> »nd p. Such a correlation might be expected if p results from inter- 
action witn the soli, w*n' (Coulomb drag and Lorentz fierce, etc,.) in which case 
-ach striae would correspond to a given charge-to-mass ratio for which the parent 
body is likely to sputter. 

Preliminary studies of the dynamics have been carried out for dust assumed 
to be released by comet Halley during its 1986 perihelion passage. The computer 
code, which also includes radiation and ion pressure is described .elsewhere 

ami MlA£Dnl-_ i _ Thft tinvn flpfaflf t rtf hll. flrtl uf rtfi { am 
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during that passage are to reinforce gravitation (counteract 8) and to deviate 
the dust towards the solar equatorial plane. For silicates corresponding to the 
cores of interstellar grains with an assumed charge as high as 10 Volt, |i w 8 ■ 0,1 
and the deviation towards the equator i3 » 20°, For "ice"-coated silicates the 
charge was estimated to 1,2 Volt, the resulting p « 8 and the deviation from the 
comet f s-orbital plane « 1.5° For larger agglomerates the interaction with the 
solar wind seera3 negligible. 

In conclusion, the dynamics of grains with high charge to mass ratio such 
as might be expected for the silicate cores resulting from the striped interstellar 
core-mantle particles, or for the bare particles, are expected to deviate sufficiently 
from central motion to distinguish them from , ‘ice' , -coated grains or agglomerates, 

8 for the parent particles are probably comparable to that for the fragments 
forming the striae in comet West. * 

We Uucn that at the limit, where the ensemble of scatterers is such that 
Mgq « f 7j~- and the scatterers are randomly located, each scatterer approaches 
its scat.ering properties as an individual particle. At this limit, 8 is the same 
for the agglomerate as for individual scatterers, 

This condition is nearly fulfilled for the "ice M ~free agglomerates especially 
V21R and TC21R, 
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5,4 Impact Data 

Analysis of lunar rock samples (Wagel et al. , 1975) reveals the existence of 
very shallow craters which are indicative of low bulk-density impacting particles, 
•p m 1 gcra“ 9 , The relationship between particle density m.d urater topology has 
futhar revealed two groups of particles stony (p * 2 - 3 gem"*) and metallic- 
density (p « 8 gem*” 3 ) meteoroids. Tho latter create deeper and much more commonly 
observed (selection effects?) craters (see also Le Sergeant and Lamy, 1978). Impact 
data from the Helios 1 spacecraft provides nn independent means of separating a 
low-density component of the interplanetary dust-complex. One of the Helios 
impact detectors was protected from direct solar radiation by a screen, The screen 
also prevented particles with density lower Chan * 1 gem” from penetrating. 

Another detector, with an open aperture pointing out of the ecliptic plane > 
registered a large excess of impacts. The data analysis revealcl a strong compo- 

*•3 

nent of low-density dust, 30 % of which may have p as low as 0,1 gem (Griin et l 
si,, 1980), Some of this dust was identified with the high-speed 8-meteoroids 
also detected by the Pioneer 8 snd 9 spacecrafts (McDonell, 1975, Berg and GrUn, 
1973), 

The higher-density "apex" particles are evidently different from the 

**•3 

agglomerates studied Here which have 1,2 > p > O.Oo gem . The "bird’s-nests" 
correspond to a low-density component orbiting lose to the invariable plane. 

Such a component world be discriminated in the Helios 1 and 2 data, A low-density 
component producing the zodiacal light would also explain the descrepancy between 
results from the Meteoroid Penetration Experiment or MPE ("beercans") experiment 
onboard Pioneer 10 and 11 spacecrafts, -showing a constant penetration rate from 
the Earth out to beyond Jupiter (Hums et al, , 1974), and the zodiacal light 
brightness observed from the same spacecrafts (for references, see Weinberg and 
Sparrow, 1978, p. ?09), On the basis of the present study it might be expected 
that the low-density 8-meteoroids are cluras of a few or single grains. The 
existence of .. higher-density component of even smaller grains (the bare grains 
in the bi-modal model for interstellar dust by Greenberg and Hong, 1973) is 
predicted. 
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5,5 Concluaiona 

It ia concluded on the basis of thiB study that the following hypothesis are 
consistent with the current observations in the interplanetary medium; 

(1) The primordial interstellar dust-complex may be described by the bi-modal 
model by Greenberg and Hong (1973), modified in the dense presolar nebula as 
described by Greenberg (1979), 

(2) Comets aggregate from the thus described interstellar dust (Whipple, 1978), 

(3) Some comets may develop an M ice M -free outer ubC (Brin and Mendis, 1979) » 

(4) The zodiacal light results from scattering of sun-light by cometary debris 
(Whipple, 1976), 

(5) 3-meteoroids might be comet and zodiacal-light particle debris. 

ihe high degree of polarization reported in cometary dust tails as far from the 

comet as the antitail, favours the "ice"~free agglomerates as representatives of 

comet debris. However, comets with large perihclie distances or "new” comets are 

possibly seeding agglomerates of primordial dust particles. These "bird's-nests" 

may be successively stripped of their volatiles as they spiral in towards the 

Sun under the action of the PoynCing-Robertson drag. The degree of polarization 

might therefore be higher for light scattered by particles closer to the Sun 

than by dust observed at large elongations, thus shifting P towards smaller 

max 

elongations. For numerous reasons (dust-dust collisions, volatile erosion, etc,.,) 
the "bird's-nests" are expected to sublimate at an increasing rate as they approach 
the Sun. The smaller or sufficiently "ice"-free debris are expelled from the 
Solar System by the radiation pressure or the solar wind and are observed as 3- 
-meteoroid8. Ocher higher-density du?* 4 , particles, such as asteroid fragments, are 
more easily observed by impact detectors and through cratering of surfaces exposed 
to the interplanetary medium, but do probably not significantly contribute to the 
zodiacal light (see also Whipple, 1976), Although the "classical" component of 
interstellar grains are thought to stream through the inner Solar System during 
part of the solar i_ycle (Gustafson and Misconi, 1979) aud might be trapped in 
bound orbits under some conditions (Gustafson, unpublished calculations) they are 
not likely candidates for the zodiacal-light particles. 

■\ 

\\ 

\\ 
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5.6 Suggestions for Further Research 

Out of the topics that suggest themselves on the basis of the present study, 
only the ones that appear to be of the most fundamental character or most 
urgent are mentioned below. No reference is made to their relative importance, 

—— Observational studies of scattering functions for comet dust 1 , 

— Investigation of the low density dust component or components detected 
by the HelioB space-probe (Griin et al., 1980). 

Renewed studies of interplanetary dust dynamics, and of the evolution 

of the zodiacal dust cloud, using the "bird's-nest" particle model 2 . 

— — Determination of the domain of validity of the suggested scattering 
approximations. 

— — Adaptation of the Finson-Probstein (1968) method to the "bird's-nest" 
type particles, for prediction of shape and structure of com et tails 2 . 


1) To the knowledge of -the author, polarimetry P(0)‘ covering an appreciable 
interval of scattering angle (0) has only been performed for comet Ikeya- 
Seki 1965 VIII. Whereas no determinations of I\.^(0) are known. Deduction of 
X (0) generally involves assumptions regarding 1- the particle number density 
distribution. However, an instrument probing the scattered intensity per unit 
volume and therefore I (0) has been suggested to fly on the forthcoming 
probe to Halley's comet? 

2) The values of the radiation pressure efficiency and the electric charge 
of the particles are of fundamental importance in evaluation of dust dynamics 
and yet often subject to wide extrapolations and guesses. 
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APPENDIX A 
— — — — : 

Scattering Relations for Arbitrary Targets and their Mirror Particles . 

a,: 

In accordance with the reciprocity theorem for vector waves, we can invert 

the time without changing the scattering problem, thus replacing H with 

-k and k Q with -k^. It follows from simple trano format ion of coordinates 

that I 2 * ■ I , I** - I ,I 2> - I and I 1 ' ■ I where I 2 * are scattered 
XX XX 22 22 12 21 21 12 

intensities from the reciprocal position. The reciprocal position is reached 

by rotating the scatterer 180° around the bisectrix of the complement to 

scattering angle. 

Assuming that for every particle there is one particle in the reciprocal 

position, I I * I T over random orientations at any given scattering angle. 
. 12 21 

bi 

Mirroring with respect to the scattering plane only changes the phase of the 
cross-polarization components, as the transformation only involves a change 
of sign of one coordinate. The attention of the reader is drawn to the fact 
that we refer to a second particle. 

a: 

A fourth related position is obtained by successively applying the two 
previous transformations a and b . 

Use of the Symmetry Relations in the Laboratory . 

Arbitrary particles . 

Only the relation described in it is valid. But with the present orientation 
mechanism the .reciprocal position can only be reached when the tilt axis 
coincides witjh the bisectrix to the complement of the scattering angle. 

Particles with an Axis of Symmetry . 

The particle appropriately reorientated is its own mirror particle. , 

The relations a, b and a all prevail. Let a vector parallel to tho sxIb 
of symmetry be our reference. Relation a implies that all scattering 
information, can be deduced as the vector sweeps the hemisphere on either 
side of the bisectrix plane, defined as the plane perpendicular to the 
scattering plane and containing the bisectrix of the complement to 8. 

Relation b implies that it is sufficient to let the vector sweep the 
hemisphere on any side of the scattering plane. Thus we need only to 
let the vector sweep a solid angle tt bound by the scattering plane and 
the bisectrix plane. Orientations in the opposite quarter of the sphere 
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are related through a, on the opposite side of the scattering plane but 
same aide of the bisectrix plane are related through b and in the fourth 
quarter through a. It is thus sufficient that measurements be made only 
in a 1/4 of the sphere. 

Particles With an Axis of Symmetry and a Plane of Symmetry Perpendicular 
to this Axis . 

In addition to a, b and a, mirroring with respect to the plane of symmetry gives 
an equivalent position. Thus, orientations in the solid angle n/2 bound by the 
scattering plane, the bisectrix plane, and a plane perpendicular to both and 
containing the geometrical center of the target, are equivalent to orien- 
tations in the opposite n/2 solid angle. By successively applying relation c* 
it is realized that it is sufficient to investigate orientations in the tt/2 
solid angle bound by the scattering plane and the bisectrix plane on -either 
side of the mutually perpendicular plane through the target. Thus for such 
particles, only. 1/8 of, the orientations fbr arbitrary particles are independent. 
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APPENDIX B 


Periodicity of the Oscillations in the Scattering Efficiency of "Bird's-nests" 
aa they are Rotated 

’ The optical path difference for light scattered by A and B (see Fig. B.1) is 
Aa x i fl varied the optical path difference varies. When; 

- ‘W 1 

-where A ia the wavelength of the incident radiation and m ifi given by eq. 3,1, 
the interference pattern has gone through n periods of oscillations. 



The period Ax is therefore defined by; 

Re{ m } b\ cos(x + Ax) - cos(x + Ax + 9) - cos(x) + cos(x + 8) J ■ A 

Note that the distances p, q and b are all in the scattering plane, b is the 
projection of the distance between two scatterers in this plane. Any pair of 
scatterers separated by the distance a and with the line joining the two centers 
making an angle a with the scattering plane Q> * cos a ) , will interfere such 
that they pass through one period of oscillation as the "bird' s-nest" is rotated 
the angle Ax from some orientation x- 

Consider a particle at the projected distance v in the scattering plane from 
the geometrical centre of the "bird' s-nest", all particles at the projected 
distance b from this particle will be on the surface of the circular cylinder of 



Fig. B.2 For each pair of interfering particles located each ot one 
of the two cylindrical surfaces theire is an orientation of 
the spherical ensemble such as the resulting interference 
pattern goes through one period of oscillation as the 
conglomerate is rotated by an angle Ax. 

radius b with the synmetry axis passing through the particle perpendicular to the 
scattering plane. As this cylinder is swept around the "b ;rd ' s-nest" (always with 
distance v between the "bird's-nest" centre and the cylinder axis), an other 
cylinder is generated also with itjt symmetry axis perpendicular to the scattering 
plane (see fig. B.2). The total number of particle pairs with the projected 

separation b is the number of particles on the first cylindrical surface times 
the number of particles on the second surface. 

A computer program generating Figs. ,8,3 looks for pairs whose individual 
projected distance are bid and where the first particle is at distance p to r 
- Ar. is therefore given by the volume V t corresponding to the first 

cylinder, times V^ t corresponding to the second cylinder, times the number density 
of particles. V is bound by two cylindrical surfaces and two spherical surfaces 
and is computed using Simpsons formula. V ^ is also bound by two cylindrical and 
two spherical surfaces, but in such a way that it can easily be computed by 
addition and substraction of simple geometrical figures. It is seen from the Figs. 
B.3 that particle pairs producing short period oscillations exists only at the 
higher scattering angles. See also section 3.6.3 p. 32 and Figs. 4.1 and 4.2 on 
p. 40 - 42. 




Total number of interfering pairs with the period 
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Fig. B.3 Total number of interfering pairs that goes through one period 
of oscillation as the ensemble is rotated an angle in the 
scattering plane from some initial orientation. The sperical 
X* 10,04 ensemble of 539 cylinders T41A is simulated at some! 
scattering angles 0. 
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APPENDIX C 

w 

Scattering from a Cloud of Randomly Orientated Independent Infinite 
Circular Cylind ers * ' 

m i mm ■ w i » n w — mi« 

A computer code to calculate scattering from single infinite circular 
cylinders at arbitrary orientation, excluding incidence along the axis of 
symmetry, as outlined by Lind (1966) and by Kerker (1969) was kindly 
provided by Dr, R. T. Wang, Note also the paper by Cohen (1980). 

The scattered light propagates along conical surfaces containing the 
H 0 vector, the symmetry axis of the cone coincides with the cylinder axis. 
Thus, only cylinders whose axis of symmetry fall in the bisectrix plane 
, (the plane perpendicular to the scattering plane and containing the bisectrix 
of the complement to 0) , will contribute to the scattering in any given 
direction. The intensities are given in components tangential and perpen- 
dicular to the cone surface. 

The geometry is shown in fig. 0,1 in transformation to the coordinates 



Fig. C.1 


The geometry for scattering by an infinite cylinder. 
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0, ij>, x in the laboratory cystem (fig, 2.1) the following relationships prevail. 

sin(0/2) “ cos(<p)sin (a/2) 

cos <U.> - s iZ/H 

sw cos 

X ■ 0/2 

The intensities contributed by individual cylinders to the total scattered 
brightness are added, as the assumption of randomness in location implies 
randomness in relative phase. 

To transform the second of Stokes parameters to a common coordinate 
system the phase relation between the four intensity components mus't be 
known as they interfere. The angular dependence of this relation is not 
yet derived, so we will only compute the total brightness. 
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Appendix D 


The Infinite Cylinder Approximation for Scattering in the Far-field b y 
Finite Cylinders 

-ft 

From Figs. 4.14 and 4.7 it seems that the angular dependence of scattering by 
finite cylinders at perpendicular incidence may be qualitatively approximated 
by that of infinite cylinders. Quantitatively the scattering is quite different. 
That off an infinite cylinder is concentrated precisely in the scattering cone 
so that the intensity at the far-field (t*» A) is inversely proportional to 
the distance r from the scstterer, while thte scattering off- finite yyliddersiis 
spread by diffraction (Lind, 1966, Greenberg, 1974) the intensity is proportional 
to v wi when *• » A. 



Fig. D.1 Experimental intensities versus scattering angle for radiation 
scattered off x m 1.88 lucite cylinders with their ‘axes of 
symmetry pointing towards 0-45° (x"45°, ^"90 ). Dash-dot is^ 
for the J, .-component of the 210:1 -elongation cylit ier extending 
* well outside the beam, daabtid curve is for the I -component 
from the same cylinder. Dotted curve is for the J ^-component 
scattered of the 4:1 -elongation cylinder C41, solid line is 
for the J 22 -component. The background radiation (appendix E) was 

not compensated for. 


i 
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The condition for scattering off an infinite cylinder towards a given 
direction is therefore that the cylinder symmetry axis is contained in the 
bisectrix plane (the plane perpendicular to the scattering plane and bisecting 
the complement to the scattering angle). Therefore, an infinite cylinder must 
be orientated perpendicular to the scattering plane in order to scatter in all 
directions 6 in the scattering plane. An infinite cylinder with the orientation 

X, »|», (see appendix C) scatter in the forward direction 0*0 and 0"2x* In 

* 

Fig, D»1 the scattering of a finite but long cylinder (extending outside the 
main lobe of both transmitting and receiving antenae) is shown. The cylinder 
is orientated such as the axis of symmetry is in the scattering plane and makes 
an angle x"44 # 6 0 with the direction of propagation of the incident radiation. 

The receiving antena was swept from 0 « 30° to 170°, as expected we see a sharp 
peak in both .uttered components end eround 8 - 2* - 89.2°. The radiation 
scattered off the finite cylinder is detectable in the whole interval of (\ short- 
ward of 0"45° part of the forward scattering peak is seen. In Fig* D.2lthe 
scattered brightness from the 4:1 -elongation cylinder averaged over "uniform" '■ 
distribution of orientations is compared qualitatively to the computed brightness 
scattered by uniformly orientated infinite cylinders, normalized at 0 * 90°. It 
is seen that the forward scattering peak extends out in side scattering thus 
seriously modifying the scattering pattern. 



{ 

i 





CD 

O 0 


Fig. 1.2 Total scattered intensities from x«1.88 lucite cylinders, veaus 

scattering angle exept 6 ■ 0°. Dots are experimental results averaged 
over "uniform" distribution of orientations of the C41 cylinder (see 
section 2.7.1). Solid curve is for the uniform distribution of ori- 
entations for infinite cylinders, computed as described in appendix 
C. 
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Appendix E 

The Background radiation in the Laboratory afc S.A.L., and' Detector Response 

Fig. B.1 shows the background radiation versus scattering angle in the laboratory 
version A, in units of the potential over the detector as measured by the lock-in 
amplifier. If not mentioned otherwise, the background radiation was suppressed as 
described in section 2.5 (p. 9- 19) to < 0.1 pV. The scattering intensity 1 in 
absolute units corresonds to *»5pV (the apparatus was calibrated using a standard 
target before and after each set of measurements and every 6 « 10 hours during the 
course of measurements) . 
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Fig. B.1 The background radiation versus scattering angle in the 
laboratory version A, when the background compensation 
waveguide was disconnected. 


In the data reduction^ a linear relationship was assumed to prevail between the 
detector unit response and the intensity of scattered radiation. This assumption 
is motivated by the measured linearity (Fig. E.2) in the range up to lOmV used 
in the measurements. 




t 
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incident power 


Fig. E.2 The response of the detector unit as a function of 
incident power. 
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